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Abstract 

"3  This  thesis  details  a  new  method  of  scheduling  project  activities  in 
project-oriented  production  systems.  The  scheduling  problem  in  a  project- 
oriented  production  system  can  be  characterized  as  a  part  of  the  hierarchical 
planning  and  scheduling  process  in  which  an  Aggregate  planning- level  allocates 
resources  among  projects  of  the  production  system  and  specifies  major  due 
dates  of  project  executions.  The  goal  of  project  scheduling  is  to  minimize  the 
project  completion  time  subject  to  resource  capacities  and  due  dates  which  are 
specified  by  the  aggregate  planning  level  Resource  capacities  for  a  particular 
project  are  time-varying  but  inflexible  at  any  particular  time.  The  method  of 
scheduling  aims  to  reflect  the  characteristics  of  work  in  ship  overhauls  in  a 
shipyard  which  is  believed  to  be  a  typical  example  of  a  project-oriented  produc¬ 
tion  system  The  method  is  designed  to  be  computationally  feasible  and 
managerially  practical. 

Although  there  are  various  methods  and  techniques  for  project  scheduling 
which  have  been  proposed  in  the  literature,  none  fulfills  all  the  needs  posed  by 
the  shipyard's  scheduling  problem.  />rV| 
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An  important  characteristic  of  work  in  a  shipyard  is  that  activities  of  pro- 
jects  may  have  variable  rates  of  resource  applications.  The  actual  duration  of 
an  activity  depends  on  the  rate  of  resource  applications.  A  model  with  variable 
intensities  for  activities  is  introduced  to  cope  with  this  work  characteristic. 
Using  the  variable  intensity  model,  heuristic  intensity  assignment  algorithms 
are  developed  for  two  different  management  policies  concerning  the  flexibility 
of  resource  redeployment.  The  intensity  assignment  algorithms  work  as  deci¬ 
sion  makers  which  decide  activity  start  times  and  intensity  levels,  ie.  resource 
allocations,  for  activities.  The  validity  of  the  algorithms  is  established  by  com¬ 
paring  their  performances  to  the  performances  of  two  methods  adapted  from 
the  methods  in  literature  in  comparative  experiments  with  arbitrarily  gen¬ 
erated  project  parameters. 

The  other  important  characteristic  of  work  in  ship  overhauling  is  the 
uncertainty  of  work  content  of  the  activities  and  also  the  uncertainty  in  the 
project  network  due  to  the  prospect  of  repair  rework  which  may  be  needed 
after  testing  of  subsystems.  To  handle  these  uncertainties,  a  simulation  model 
which  reflects  the  stochastic  nature  of  work  and  uses  the  intensity  assignment 
algorithms  as  a  decision  maker  is  introduced. 

Finally,  the  intensity  assignment  algorithms  and  the  simulatidn  model  are 


used  to  generate  several  schedules  for  an  actual  large  ship  overhaul  project  as 
an  implementation  of  the  method  proposed.  Results  are  discussed. 
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1.  INTRODUCTION 


Scheduling  large  projects  is  a  very  important  task  in  project- oriented  pro¬ 
duction  systems,  such  as  shipyards  and  airplane  construction  firms.  We  will 
refer  to  a  production  system  that  completes  a  number  of  different  types  of 
large  projects  on  a  repetitive  basis  as  their  outputs  as  a  project-oriented  pro¬ 
duction  system.  Such  systems  are  characterized  by  inflexible  capacities  for 
scarce  resources,  particularly  various  types  of  skilled  trade  labor.  These 
resources  must  be  utilized  efficiently  in  multiple,  concurrent  projects  in  order 
to  maximize  throughput  of  the  system. 

The  scheduling  of  individual  projects  has  been  given  considerable  attention 
from  both  theoretical  and  practical  people  because  of  its  substantial  financial 
impact.  In  general,  the  cost  of  large  projects  is  very  high  and  even  a  small 
reduction  in  project  time  may  realize  substantial  savings  for  a  company. 

A  number  of  methods  have  been  developed  and  a  large  amount  of  computer 
software  has  been  marketed  in  the  last  three  decades  for  the  project  schedul¬ 
ing  problem.  In  general,  methods  developed  are  tailored  to  schedule  the  pro¬ 
jects  for  construction  firms  where  the  resource  capacities  are  flexible  and  the 
projects  undertaken  by  those  companies  are  not  similar  to  one  another. 

The  scheduling  problem  in  a  project-oriented  production  system  can  be 
characterized  as  a  part  of  the  hierarchical  planning  and  scheduling  process  in 
which  an  aggregate  planning  level  allocates  resources  among  projects  of  the 
production  system  and  specifies  major  due  dates  of  project  executions.  The  goal 
of  project  scheduling  is  to  minimize  the  project  completion  time  subject  to  the 
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resource  allocations  and  due  dates  which  are  specified  by  the  aggregate  plan¬ 
ning  level.  The  resource  allocations  to  a  particular  project  serve  as  resource 
capacities  for  the  project.  These  capacities  are  time-varying  but  inflexible  in  a 
particular  time. 

The  lack  of  a  scheduling  method  that  captures  the  needs  of  project- 
oriented  production  systems  led  to  the  development  of  the  method  described 
herein  for  scheduling  projects  in  such  systems.  Chapter  2  describes  the 
scheduling  environment  in  project-oriented  production  >>stems.  Differences 
between  these  systems  and  one-time  project  systems  are  identified. 

A  review  of  project  scheduling  methods  is  given  in  Chapter  3.  The  review 
reflects  the  historical  development  of  the  project  scheduling  problem.  Resource 
Unconstrained  Project  Scheduling.  Resource  Constrained  Project  Scheduling, 
and  Generalized  Project  Scheduling  methods  are  reviewed. 

Modeling  assumptions  for  a  new  planning  and  scheduling  system  in  a 
project-oriented  production  environment  are  given  in  Chapter  4.  A  generalized 
model  for  Resource  Constrained  Project  Scheduling  is  introduced  in  Chapter  5. 
In  this  model,  variable  resource  application  levels  which  we  call  variable  inten¬ 
sities  for  each  activity  are  used.  A  linear  activity  analysis  model  is  used  to 
characterize  variable  intensities  for  activities.  Intermediate  product  transfers 
between  activities  are  analyzed. 

Chapter  6  provides  the  development  of  intensity  assignment  algorithms  for 
two  different  management  policies.  The  intensity  assignment  algorithms  work  as 
decision-makers  which  decide  activity  start  times  and  intensity  levels.  In  the 
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first  policy,  management  does  not  allow  reduction  in  the  level  of  resources 
which  is  committed  to  an  activity  once  it  is  started,  although  it  allows  an 
increase  in  this  level  while  the  activity  is  in  progress  (Upgrading-Only  Algo¬ 
rithm).  In  the  second  policy,  management  allows  reductions  in  the  level  of 
resources  that  is  committed  to  an  activity  as  well  as  increases  in  it  (Upgrading- 
Downgrading  Algorithm).  Both  algorithms  use  the  same  procedure  to  handle 
time  varying  resource  capacities  which  is  explained  in  section  6.S. 

Chapter  7  demonstrates  the  validity  of  the  intensity  assignment  algo¬ 
rithms.  Two  other  scheduling  methods  from  the  literature  and  the  intensity 
assignment  algorithms  are  all  used  to  schedule  arbitrarily  generated  projects. 
(The  two  methods  used  in  comparison  are  described.)  The  procedure  for  arbi¬ 
trary  network  generation  is  summarized.  Comparative  results  of  experiments 
are  analyzed. 

Chapter  8  proposes  new  methods  for  handling  stochastic  elements  in  pro¬ 
ject  scheduling,  focusing  on  stochastic  resource  requirements  and  network 
structure.  Stochastic  parameters  are  described.  Treatments  for  these  parame¬ 
ters  during  simulations  are  explained. 

Finally,  intensity  assignment  algorithms  are  utilized  to  simulate  a  ship 
overhaul  project  and  results  are  presented.  The  implications  of  the  results  for 
shipyard  management  are  given  in  Chapter  9. 


2.  SCHEDULING  IN  PROJECT-ORIENTED  PRODUCTION  SYSTEHS 


Project-oriented  production  systems  produce  large  similar  projects  as  its 
products.  In  order  to  accomplish  the  completion  of  one  unit  of  such  a  product, 
a  large  number  of  interrelated  activities  are  performed  in  sequence.  Each 
activity  requires  work  from  various  trade  labor  types.  As  an  example  of  such  a 
production  system,  a  shipyard  for  overhauling  naval  ships  can  be  given.  At  any 
given  time  ten  to  fifteen  ships  may  be  in  various  of  stages  of  overhaul  in  such  a 
shipyard.  Overhaul  of  different  ships  in  the  same  class  have  very  similar  tasks. 
Overhaul  of  a  ship  may  consist  of  thousands  of  activities.  Each  activity  may 
require  work  effort  from  hundreds  of  different  trade  labor  shops. 

2.1.  Production  System  versus  One-time  Projects 

We  view  a  production  system  as  a  system  which  uses  a  stream  of  inputs 
(labor  of  different  types  .  raw  material,  money,  production  facilities  including 
equipment  and  machinery)  in  order  to  produce  some  design  outputs  (projects 
in  our  case).  In  a  production  system  there  are  permanent  facilities  and  per¬ 
manent  resources  (trade  labor)  in  the  system.  Those  resources  are  utilized  by 
the  projects  which  are  in  progress.  Note  that,  even  if  there  is  no  production, 
the  cost  of  those  resources  and  facilities  must  be  absorbed  by  the  system. 

The  project  cost  for  the  production  system  largely  depends  on  project  time 
and  allocation  of  the  system  costs  among  projects.  The  objective  of  planning 
and  scheduling  the  projects  and  project  activities  in  such  a  system  is  to  take  on 
and  finish  the  projects  as  early  as  possible  within  the  given  resource  capacities. 
Equivalently,  the  objective  is  to  maximize  the  throughput,  ie,  the  number  of 
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projects  completed  by  the  system  within  a  given  (rolling)  time  horizon. 

We  will  refer  to  projects  that  are  executed  once  (or  perhaps  a  small 
number  of  times)  as  one-time  projects.  In  systems  performing  one-time  projects 
the  facilities  and  resources  which  will  be  used  in  a  project  are  obtained  tem¬ 
porarily.  The  acquisition  or  rental  of  facilities  and  resources  are  made  after  a 
decision  is  made  to  take  on  the  project.  Timing  and  amounts  of  resources  that 
will  be  acquired  for  the  project  can  be  planned  in  advance.  Also,  supplementary 
resources  may  be  acquired  during  the  project  execution  if  any  bottleneck  may 
arise.  Most  of  the  construction  projects  (building  bridges,  dams,  etc.)  can  be 
considered  as  one-time  projects  since  construction  companies  do  not  hold  sub¬ 
stantial  amounts  of  permanent  resources.  For  one-time  projects,  project  cost 
largely  depends  on  the  cost  of  the  resources  that  are  temporarily  acquired.  The 
objective  of  planning  and  scheduling  of  one-time  projects  is  to  minimize  the 
cost  of  the  project  while  keeping  the  promised  project  due  dates. 

The  research  reported  here  is  the  development  of  some  techniques  to 
assist  project  scheduling  in  project-oriented  production  systems.  The  research 
is  motivated  by  a  study  of  ship  overhaul  scheduling  within  a  naval  shipyard. 

2.2.  Single  Project  versus  Multi-Project  Environments 

In  a  single  project  environment,  the  production  system  executes  only  one 
project  at  a  time.  Total  resources  of  the  production  system  are  available  to 
activities  of  the  project  which  is  in  execution. 

In  a  multi-project  environment,  there  are  simultaneous  or  overlapping  pro¬ 
jects  which  are  competing  for  the  same  resources.  For  scheduling  purposes. 


one  might  consider  treating  the  collection  of  projects  as  one  big  project.  This 
can  be  practical  if  the  projects  are  not  too  big;  otherwise  the  scheduling  prob¬ 
lem  becomes  intractable. 

Boysen  (1982)  and  Leachman  and  Boysen  (1984)  have  proposed  manage¬ 
ment  in  large  multi-project  environments  using  multiple  levels  of  planning  and 
scheduling  effort.  Planning  and  scheduling  in  a  multi-project  production 
environment  can  be  viewed  from  the  perspective  of  both  project  management 
and  resource  management.  Effective  planning  and  scheduling  will  only  be 
achieved  if  those  two  perspectives  are  properly  linked  together. 

Planning  projects  and  resources  is  a  complex  task  in  multi-project  produc¬ 
tion  environments,  in  which  hierarchical  levels  of  planning  must  necessarily 
exist.  Boysen  (1982)  proposes  three  levels  of  planning  in  the  planning  process 
in  such  a  system. 

These  levels  are  as  follows: 

1.  Inter-Project  Planning 

2.  Individual  Project  Planning 

3.  Project  Scheduling 

Coing  from  higher  levels  to  lower  levels,  additional  constraints  are  imposed, 
the  planning  horizon  becomes  shorter,  there  is  more  specific  information,  and 
the  degree  of  uncertainity  decreases.  Explanations  of  decisions  involved  in  each 
planning  level  are  provided  in  section  2.3. 
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2.3.  Projecl  Planning 


Inter-Project  Planning 

Long  range  planning  is  needed  to  fit  new  projects  into  the  system.  The  tim¬ 
ing  and  compression  of  the  projects  largely  depends  on  the  resource  loadings  of 
each  project.  A  new  project  “fits  in"  if  a  consistent  and  feasible  plan  for  man¬ 
ning  the  new  project  and  all  other  projects  is  developed.  Expected  due  dates  of 
some  important  milestones  of  the  project  must  be  estimated.  On  the  basis  of  the 
inter-project  plan,  time-dependent  project  progress  workload  forecasts  are 
given  to  resource  managers. 


Individual  Project  Planning 

The  aim  of  this  level  is  to  track  and  control  project  execution  to  keep  pro¬ 
jects  on  schedule.  Projects  are  described  in  more  detail  in  this  level.  Certain 
project  due  dates  will  have  been  established  at  the  multi  project  planning  level. 
Additional  milestones  are  then  defined  to  indicate  start  and  finish  of  important 
phases  of  the  project. 

A  detailed  network  of  activities  is  developed  to  observe  the  progress 
between  milestones  more  accurately.  Duration  and  resource  requirements  for 
each  activity  are  estimated.  Using  these  estimates  for  durations,  predictions  of 
milestone  actualization  dates  can  be  made  more  accurately.  Flexibility  in  per¬ 
forming  the  project  activities  can  be  determined. 

Total  resource  requirements  of  the  project  activities  over  the  project  dura¬ 
tion  are  forecasted  more  accurately.  The  feasibility  of  the  allocation  of 


resources  to  the  project  in  the  inter*project  analysis  may  be  verified  and/or 
modifications  of  input  parameters  for  inter-project  analysis  may  be  made.  Using 
detailed  activity  networks  it  is  easier  and  more  accurate  to  track  and  control 
the  project. 

Project  Scheduling 

This  level  involves  day-to-day  operating  decisions.  The  detailed  activity- 
network  is  used  to  establish  the  priorities  among  the  project  activities.  The 
resources  preallocated  to  the  project  by  higher-level  planning  are  allocated 
among  the  activities  of  the  project  using  the  scheduling  priorities.  Scheduling 
of  activities  is  constrained  by  their  precedence  relations,  preassigned  due  dates 
and  preallocation  of  resources  of  each  type.  A  daily  schedule  shows  the  activi¬ 
ties  that  will  be  in  execution  and  the  amounts  of  resources  used  by  them. 

The  project's  execution  status  is  monitored  and  updated.  The  status  can  be 
measured  in  terms  of  progress  in  activities.  The  initial  durations  and  resource 
requirements  of  activities  are  revised  if  necessary.  Delay  in  activities  is 
identified. 

The  research  reported  here  deals  with  the  project  scheduling  level. 
Improvement  in  the  scheduling  level  will  facilitate  improvement  of  the  planning 
process  in  all  levels.  It  serves  to  complement  the  work  of  Boysen  (1982)  and 
Leachman  and  Boysen  (1984)  on  the  methodology  for  inter-project  planning. 

2.4.  Example  of  Industry  Practice:  Shipyard's  Approach 

A  shipyard  for  overhauling  naval  ships  is  examined  as  an  example  of  a 


project-oriented  production  system. 

There  are  three  hierarchical  planning  levels  for  project  planning  and 
scheduling  in  the  shipyard. 

1.  Inter-Project  Planning 

2.  Project  Planning 

3.  Project  Control 

The  inter-project  planning  level  decides  the  timing  of  future  ship  overhauls 
for  a  horizon  of  a  few  years.  The  decision  process  can  be  described  as  follows. 

The  shipyard's  planning  stall  has  rough  information  of  work  requirements 
of  future  overhauls.  The  source  of  the  information  is  primarily  past  experience, 
since  ships  overhauled  belong  to  specific  classes  and  ships  in  the  same  class  are 
very  similar.  Each  ship  overhaul  has  several  alternative  loading  curves  for  total 
manpower.  The  shapes  of  those  curves  are  examples  from  past  history  (when 
available),  scaled  to  the  requirements  of  the  upcoming  overhaul.  The  staff  ini¬ 
tially  selects  a  loading  curve  for  each  future  overhaul.  Using  these  curves,  a 
rough-cut  capacity  analysis  is  performed  to  establish  "reasonable"  start  and 
due  dates  for  future  overhauls.  The  loading  curves  (ie.,  load  histories  of  the 
sum  of  all  trade  labor  shops)  for  ongoing  and  future  overhauls  are  juxtaposed 
and  tabulated  according  to  the  trial  start  and  finish  dates  for  each  overhaul. 
The  resulting  total  load  is  compared  to  total  shipyard  manpower.  Ideally,  the 
selection  of  loading  curves  would  be  made  so  that  the  sum  of  manpower  load 
across  overhauls  fits  the  "capacity"  of  the  yard  within  "reasonable"  tolerance. 
Where  major  discrepancies  exist,  curves  are  shifted  or  substituted  until  one 


comes  up  with  "reasonable"  start  and  due  dates  for  future  projects. 

Once  these  rough  start  and  finish  times  of  projects  are  decided  by  the 
planning  staff,  the  workload  forecasting  department  prepares  "Workload  Fore* 
casts"  corresponding  to  the  total  manpower  curve  selected  for  each  overhaul. 
Summing  up  each  resource  across  overhauls  according  to  the  selected  resource 
loading  patterns  gives  the  workload  forecast  for  that  resource  for  that  planning 
horizon. 

Forecasts  of  total  resource  workload  and  the  resource  workload  of  each 
overhaul  are  given  to  resource  managers  so  that  they  can  make  adjustments  of 
their  work*force.  Ordinarily,  no  major  adjustments  in  work-force  levels  are  rou¬ 
tinely  allowed  except  the  timing  of  labor  vacations. 

As  the  starting  time  for  an  upcoming  overhaul  draws  nearer,  project  plan¬ 
ning  is  undertaken,  whereby  the  work  requirements  are  defined  and  estimated 
in  more  detail.  A  critical  path  network  is  developed  and  activity  durations  are 
estimated.  Due  dates  for  major  project  events  are  established  by  judgement  of 
senior  management. 

The  Project  Control  level  aims  to  track  and  control  the  progress  of 
overhauls  in  order  to  finish  the  overhauls  at  planned  due  dates,  and  within 
resource  capacities. 

Each  overhaul  is  assigned  a  project  manager  who  is  responsible  to  track 
and  control  the  project.  The  scheduling  department  helps  the  project  manager 
to  control  the  project.  It  prepares  a  detailed  activity  network  for  the  project.  It 
does  ordinary  critical  path  method  calculations  for  the  overhaul  by  using  the 
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due  dates  for  major  events  specified  in  project  planning  as  constraints  to 
schedule  the  activities.  Schedules  developed  here  are  not  necessarily  con¬ 
sistent  with  resource  capacities,  since  CPM  calculations  are  made  assuming 
that  there  are  enough  resources  in  the  system.  The  scheduling  department 
updates  the  CPM  calculations  periodically. 

For  each  activity,  the  shop  with  the  most  work  required  is  designated  as 
the  lead  shop  for  the  activity.  Each  trade  labor  shop  has  a  manager  who  assigns 
labor  among  activities  of  projects.  Every  shop  manager  is  furnished  with  a  list 
of  priorities  among  the  activities  on  the  basis  of  least  total  float  of  activities 
from  CPM  calculations.  The  shop  manager  allocates  manpower  from  all  shops 
among  available  activities  for  which  his  shop  is  designated  as  lead  shop.  In  gen¬ 
eral.  activities  can  be  manned  in  many  different  ways.  The  lead  shop  managers 
are  the  ones  who  decide  the  way  of  manning  the  activities. 

The  project  manager  is  responsible  for  controlling  the  project  execution.  If 
any  critical  activity  is  delayed  because  a  resource  needed  by  the  activity  is  not 
allocated  by  a  shop  manager,  the  project  manager  reports  it  to  top  manage¬ 
ment.  Top  management  may  resolve  the  conflict  between  project  managers  and 
shop  managers  by  establishing  priorities  between  projects  which  are  competing 
for  the  resource.  It  may  overrule  the  allocations  made  by  the  shop  manager  in 
the  case  that  the  project  manager  successfully  appeals  to  top  management. 


2.5.  Special  Characteristics  of  York 


Flexible  Loading  of  Resources 

Many  activities  of  projects  in  tbe  shipyard  may  be  manned  in  many 
dilTerent  ways.  Lead  shop  managers  have  a  large  amount  of  flexibility  in  man¬ 
ning  activities.  The  total  resources  required  by  the  activity  must  be  satisfied 
but  the  rate  and  pattern  of  resource  application  is  decided  by  the  shop 
manager.  He  has  tbe  flexibility  of  allocating  as  many  men  as  is  practical  if  he 
thinks  the  job  is  urgent  or  he  may  allocate  a  smaller  crew  otherwise.  The  actual 
duration  of  the  job  depends  on  the  rate  of  resource  applications. 

Uncertain  York  Content 

Each  activity  of  a  project  has  authorized  work  from  the  shops  whose  ser¬ 
vices  are  required.  The  authorized  work  required  by  an  activity  from  the  shops 
is  defined  in  terms  of  tbe  tasks  to  be  accomplished,  but  in  actuality  there  is 
considerable  uncertainty  in  tbe  total  aunount  of  resources  required,  ie,  the 
total  man-hours  to  accomplish  the  tasks  required  by  the  job.  Historical  data  on 
total  man-hours  from  various  shops  required  by  each  activity  is  available  from 
which  distributions  for  total  man  hours  can  be  developed. 

Besides  authorized  work  from  the  shops  for  tbe  project  activities,  the  need 
for  rework  or  previously  unauthorized  work  may  be  discovered  during  testing  of 
the  subsystems  that  have  been  worked  on  before.  Testing  of  the  subsystems  is  a 
substantial  part  of  the  overhauling  project.  In  the  shipyards  rework  of  author¬ 
ized  work  or  newly-discovered  (unauthorized)  work  occurs  frequently. 


rV 


vrvtui.ru 


Needs  of  Shipyard  Management 

Shipyard  management  needs  a  system  design  for  project  planning  and 
scheduling  which  reflects  the  characteristics  of  work  in  overhauling  which  are 
discussed  above  and  which  fits  in  with  inter-project  planning.  Proposed 
methods  should  be  computationally  feasible  and  managerially  practical.  The 
study  reported  here  is  motivated  by  the  needs  of  shipyard  management.  The 
discussion  of  the  proposed  system  is  deferred  until  Chapter  4  in  order  to  review 
previously-developed  scheduling  methods  (Chapter  3). 


3.  REVIEW  OF  PROJECT  SCSDULMC 


3.1.  Ovtntev 

In  general,  a  large  complex  project  mvohrvs  a  large  number  of  component 
activities  that  are  dependent  on  each  ether  by  specified  precedence  require* 
ments.  Scheduling  the  component  activities  in  order  to  achieve  a  certain  goal 
is  the  main  task  of  project  management.  The  goal  of  scheduling  may  be.  com* 
pleting  the  project  by  a  specific  deadline,  or  minimising  the  cost  of  the  project 
in  order  to  meet  the  specified  target  due  date,  or  minimising  the  project  com* 
pletion  time.  etc.  The  constraints  imposed  on  scheduling  of  a  project  may  also 
vary  depending  on  the  specific  requirements  of  the  project. 

Methods  of  scheduling  may  vary  according  to  their  specific  goal  and  con* 
straints  imposed  on  them  but  they  have  an  important  feature  in  common:  all  of 
them  are  based  on  the  representation  of  project  as  a  network  of  activities. 
Therefore,  they  often  referred  to  as  network  analysis,  network  planning,  and 
network  scheduling  The  most  widely  used  names  are  CPU  and  PERT,  which 
stand  for  two  techniques  that  evolved  in  the  late  fifties,  and  they  will  be  covered 
in  section  3.2. 

Activity  Networks 

A  project  can  be  represented  by  an  activity  network  in  which  nodes 
represent  activities  and  arcs  represent  precedence  relations  between  activities. 
This  representation  is  called  as  Activity-on-Node.  Activity-on-Arc  representa¬ 
tion  is  also  used  but  it  requires  introduction  of  dummy  activities  to  show  some 
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precedence  relations  which  is  not  needed  A-O-N  representations'!- 

Projects  are  decomposed  into  activities  using  several  criteria.  Resource 
use,  work  location,  and  technical  system  are  common  criteria  for  project 
decomposition. 

An  activity  is  a  part  of  a  project  which  has  beginning  and  ending  points, 
clearly  identifiable,  which  are  called  events.  It  has  a  duration,  which  is  the  time 
interval  elapsed  between  the  moment  the  activity  starts  and  the  moment  the 
activity  is  completed.  It  consumes  resources  of  different  types.  It  is  allowed  to 
start  any  time  after  its  predecessors  are  completed. 

Review  of  project  scheduling  methods  for  the  resource-unconstrained  case 
will  be  presented  in  section  3.2.  Even  though  the  assumption  that  resources 
are  unconstrained  is  not  a  reasonable  one,  the  concepts  developed  for  this  case 
are  largely  used  in  the  resource  constrained  project  scheduling  methods,  which 
are  reviewed  in  section  3.3.  Section  3.5  provides  a  review  of  generalized  project 
scheduling  methods.  Emphasis  will  be  given  to  the  methods  that  generalize  the 
scheduling  problem  in  terms  of  resource  application  flexibility  and  stochastic 
elements  added  to  the  problem.  Section  3.5.1  and  3.5.2  covers  these  two  classes 
respectively. 

3.2.  Resource  Unconstrained  Project  Scheduling  Methods 

There  are  two  well  known  methods  for  planning  and  scheduling  the  project 
for  the  case  that  resources  are  unconstrained.  These  are  CPU  (  Critical  Path 

t  Other  advantages  for  use  of  A-O-N  representation  are  given  by  Hoder  and  Philips()964) 
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Method  )  and  PERT  (  Program  Evaluation  and  Review  Technique).  They  both  use 
similar  network  representations  and  activity  scheduling.  Since  the  two 
methods  were  developed  independently  ,  they  bore  the  mark  of  environments 
for  which  they  were  created. 

CPU  was  a  result  of  a  joint  effort  of  DuPont  Company  and  Remington  Rand 
Univac  Division  originally  aimed  at  better  planning  and  control  of  the  overhaul 
and  maintenance  of  chemical  plants.  The  project  was  started  in  1957  and  the 
CPM  method  was  released  in  1959  (  Kelley  and  Walker  ). 

PERT  was  developed  1957*1956  by  a  research  team  set  up  by  the  Navy  Pro¬ 
ject  Office.  The  methodology  was  published  in  1959  (  Malcolm  et  al.).  PERT  was 
developed  for  and  has  been  used  most  frequently  in  aerospace  industries,  ie, 
R&D  types  of  programs  (Elmaghraby  (1977)).  These  industries  are  relatively 
new;  their  technology  is  rapidly  changing,  and  their  products  are  non-standard. 
Thus  PERT  had  to  cope  with  the  uncertainties  that  accompany  R&D  activities. 
Hence,  PERT  regards  the  total  project  duration  as  a  random  variable,  and  per¬ 
forms  probabilistic  calculations  in  order  to  characterize  it. 

CPM  was  developed  in  the  context  of  routine  operations  whose  durations 
were  more  or  less  well  established.  Therefore,  CPM  is  basically  deterministic. 
The  major  consideration  of  the  originator  of  CPM  was  to  reduce  the  total  time 
during  which  the  facility  had  to  be  shut  down  for  overhaul. 

3.2.1.  Critic*]  Path  Method 

The  Critical  Path  Method  is  a  technique  for  determining  the  ''feasible” 
schedules  (  See  Moder  &  Philips  (1970)  or  Elmaghraby  (1977).).  The  duration  of 
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each  activity  A ,,  denoted  by  d*.  is  assumed  prespecified.  It  is  also  assumed  At  is 
not  interruptible.  Let  55=(Sj . Sjv)  be  a  schedule  of  activity  start  times. 

Let  Aj  be  a  successor  of  A*.  If  Aj  may  not  start  until  At  has  finished  the  fol¬ 
lowing  inequality  must  hold: 

^  *  ft*  *  Sj 

Given  all  precedence  relations  relating  the  activity  start  times  ,  the  earliest 
•tart  time  .  £5*  .  of  each  activity  At  is  computed  by  moving  forward  through  the 
network.  The  earliest  possible  project  completion  time  can  be  determined,  ie. 
the  earliest  finish  time  of  the  last  activity  in  the  project.  Using  this  date  ,  or 
some  specified  completion  time  exceeding  it.  one  can  move  backward  through 
the  network,  thereby  computing  the  latest  start  time  of  each  Af  The  early 

schedule  is  denoted  herein  by  £5=(£S1 . £S*)  and  the  late  schedule  by 

LS=(LS1  LSk)  .  Every  feasible  schedule  5  thus  satisfies, 

ESt^StSLSt  i- 1 . N 

The  slack  of  activity  i  is  defined  as  TSt  =  £5*  -  ES^.  It  is  the  amount  of 
time  that  the  start  of  At  can  be  delayed  beyond  £5*  if  its  successors  operate 
late  and  its  predecessors  early.  The  activities  with  zero  slack  form  the  critical 
pathos)  of  the  network. 

3.2.2.  PERT  Method 

As  mentioned  earlier,  PERT  and  CPM  use  the  same  network  representation 
and  activity  scheduling  methods.  PERT  assumes  activity  durations  are  uncer¬ 
tain.  For  each  activity  in  the  project  network,  rather  than  trying  to  estimate 
directly  the  expected  performance  time  of  an  activity,  PERT  developers  chose 
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to  estimate  three  durations  instead: 

*  The  optimistic  duration,  o  ,  is  the  shortest  possible  time  to  complete  the 
activity,  assuming  most  favorable  conditions. 

*  The  pessimistic  duration,  p  ,  is  the  maximum  time  that  would  be  necessary 
to  perform  the  activity  under  most  unfavorable  conditions. 

*  The  most  likely  duration,  m. 

To  use  those  values  one  must  specify  the  kind  of  probabilistic  distribution 
associated  with  the  activity  duration.  Clark  (1962)  proposed  to  use  beta  distri¬ 
butions  for  this  purpose.  His  proposal  has  become  extensively  accepted.  With 
this  assumption,  the  estimated  average  duration.  //,  and  .  estimated  standard 
deviation  of  activity  duration,  o  are: 


Having  the  expected  activity  durations  one  can  perform  scheduling  compu¬ 
tations  in  section  3.2.1.  For  this  case,  we  are  mostly  interested  in  a  "forward 
pass"  computation  that  yields  the  expected  earliest  start  time  of  each  A*. 

Some  milestones,  ie  .  important  events,  for  the  project  may  have  pre¬ 
assigned  scheduled  occurrence  times.  Let  5%  be  the  scheduled  time  for  start  of 
A*.  Then,  knowing  £3*  and  Si,  the  question  is  to  what  extent  can  we  assert  that 
the  scheduled  time  will  met.  This  is  simply  the  problem  of  determining  a 
confidence  interval  around  estimated  length  of  a  particular  path  in  the  net¬ 
work.  A  path  is  a  concatenation  of  activities.  Therefore,  its  length,  as  the  sum  of 
uncertain  durations,  is  itself  a  random  variable.  Probability  theory  results  for 
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sums  of  random  variables  enables  us  to  estimate  the  probabilistic  characteris¬ 
tics  of  path  length,  ie  ,  path  duration  and  variance. 

In  general  PERT  tends  to  yield  optimistic  estimates  for  expected  path 
lengths  and  for  the  probabilities  of  achieving  events  by  their  schedule  dates. 
This  is  mainly  because  the  assumption  that  the  critical  path  is  fixed  and  given 
once  for  all  is  wrong  in  principle  due  to  the  variability  in  activity  durations.  An 
extensive  investigation  of  this  bias  and  of  the  assumptions  underlying  the  PERT 
calculations  has  been  made  by  Mac  Crimmon  and  Ryavec  (1961).  Other  authors 
like  Clark(196l).Fulkerson  (1962).  Clingen(1967).  Martin  (1965),  Elmaghraby 
(1967),  and  Robillard  and  Trahan(1977).  to  mention  a  few  .  have  studied  the 
problem  in  search  of  improved  estimates  of  path  length  distribution  functions. 
However,  the  dependency  among  path  lengths  renders  the  analytical  treatment 
of  correcting  for  the  PERT  bias  quite  difficult. 

Monte  Carlo  simulation  represents  an  alternative  approach  to  this  problem. 
Van  Slyke  (1963)  uses  this  simulation  procedure  to  establish  the  unbiased  dis¬ 
tribution  of  the  project  duration.  In  principle,  each  simulation  runs  consists  of 
drawing  at  random,  for  each  project  activity,  a  duration  time  from  the  associ¬ 
ated  distribution  function,  and  then  computing  the  resulting  project  schedule. 
After  a  suitable  number  of  runs,  one  can  draw  estimates  of  the  mean  and  vari¬ 
ance  of  the  project  duration. 

3.3.  Resource  Constrained  Project  Scheduling  Problem 

Consider  the  problem  of  scheduling  the  activities  within  a  project  subject 
to  resource  use  constraints.  The  literature  addresses  a  number  of  related 
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resource  constrained  scheduling  problems  (Herroelen  (1972),  Patterson  (1973). 
Davis  (1973).  Cooper  (1976).  Elmaghraby  (1977).  Stinson  (1978).  Talbot  (1982). 
Hax  and  Candea  (1984)). 

Time /Cost  Tradeoff  problem:  The  resource  cost  of  performing  each  activity 
is  expressed  as  a  function  of  its  duration.  The  problem  is  to  find  a  schedule 
which  minimizes  the  total  cost  of  the  project.  Linear  and  convex  cost-duration 
functions  are  commonly  used  in  solution  methods  (  Fulkerson  (1961).  Kelley 
(1961).  Moder  and  Phillips  (1970),  Phillips  and  Dussouky  (1977)  ).  General  more 
complicated  discontinuous  cost  functions  can  also  be  handled  for  very  small 
problems  (  Meyer  and  Shaffer  (1963),  Moder  and  Phillips  (1970)  ). 

Resource  Loading  Analysis  .For  a  given  project  schedule,  one  produces  time 
charts  regarding  the  usage  of  various  of  resources  required  by  the  project 
activities.  Charts  may  be  used  in  checking  against  resource  availabilities  to 
establish  whether  the  project  is  feasible.  Or,  an  analysis  of  alternative  loading 
profiles  resulting  from  different  project  schedules  may  suggest  ways  to  further 
improve  resource  utilization. 

Resource  Leveling :  Costs  are  assumed  to  be  a  function  of  the  peak  use  of 
each  resource.  Then,  the  problem  is,  given  a  prespecified  project  completion 
time  and  resource  usage  by  activities,  to  find  the  schedule  that  minimizes  the 
resource  capacity  costs.  In  general,  as  the  combinatorial  nature  of  the  problem 
precludes  optimal  solutions,  leveling  procedures  are  heuristically  based  (Bur¬ 
gess  and  Killebrew  (1962).  Levy  et  al.(1962).  DeWitte  (1964).  Moder  and  Phillips 
(1970).  Leachman  (1963)). 


Resource  Constrained  problems:  For  given  resource  capacities,  the  project 
duration  is  to  be  minimized.  Since  the  problem  we  try  to  deal  with  falls  into  this 
category,  special  attention  will  be  given  to  this  type  of  problem  in  the  review. 

3.4.  Solution  Methods  for  Resource  Constrained  Problem 

In  the  resource  constrained  problem  is,  one  allocates  the  limited  resources 
to  competing  activities  in  order  to  achieve  a  schedule  which  results  in  the 
minimal  project  completion  time  or  a  schedule  which  has  minimum  cost.  The 
problem  assumes  activities  have  fixed  durations  and  fixed  resource  require¬ 
ments. 

The  problem  can  be  conceptualized  by  a  mathematical  formulation 
adapted  from  Talbot  and  Patterson  (1978). 

The  following  variables  are  used  in  the  formulation: 

Ft  Scheduled  finish  time  for  activity  i;  a  project  schedule  is  defined  by  a 

complete  set  of  Ft  . 

eft  Activity  duration. 

Pi  Set  of  immediate  predecessors  of  activity  i. 

t  Time,  assumed  to  be  divided  into  integral  periods,  indexed  by  t  =  1 T. 

At  Set  of  activities  in  progress  at  period  t  . 

a*  Usage  rate  of  resource  type  k  by  activity  i. 

Amount  of  resource  type  k  available  in  time  period  f  . 

For  the  objective  of  minimizing  the  project  duration,  the  mathematical 


model  is  as  follows: 


Minimize  Max  [ft] 

(m  twtrk 

Subject  to: 

Max  [ft]  +  iL  *  ft  itnxtwerk 
ftPt  ’ 

2® . t 

tU| 

The  combinatorial  nature  of  the  problem  is  apparent  if  we  consider  4. 
Different  project  schedules  may  have  different  sets  of  4-  A  number  of  pro¬ 
cedures  are  proposed  to  solve  this  problem.  One  can  easily  classify  the  solution 
methods  into  two  types,  one  being  exact  methods  (  methods  which  produces 
optimal  schedules),  and  the  other  being  heuristic  methods. 

3.4.1.  Exact  Methods 

Integer  Programming  and  various  Enumeration  Methods  are  used  to  find  an 
optimal  solution  to  the  problem. 

Integer  Programming 

0-1  integer  variables  are  used  to  identify  the  elements  of  the  active  activity 
set  4,  yielding  an  integer  programming  model.  The  problem  is  still  hard  to 
solve.  In  general  the  size  of  the  problem  that  oan  be  handled  by  an  integer  pro¬ 
gram  is  found  to  be  very  limited  (  Pritsker  at  al.  (1969),  Wiest  and  Levy  (1977)  ). 
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Enumeration  Method* 

One  approach,  that  can  easily  be  seen,  is  to  enumerate  all  possible 
schedules,  and  select  the  best.  However  simple  enumeration  is  computationally 
infeasible.  A  number  of  people  have  tried  to  cut  the  size  of  the  enumeration 
either  by  introducing  some  heuristic  that  reaches  the  optimal  solution  faster, 
or,  exploiting  the  special  tree-like  structure  of  the  project  scheduling  problem 
which  allows  sequential  decisions  (  Mueller-Mehrback  (1967),  Johnson  (1967), 
Davis  (1969),  Davis  and  Heidorn  (1971),  Talbot  and  Patterson  (1978)  ).  Thus 
various  implicit  enumeration  methods  have  been  devised  to  find  optimal 
schedules.  These  methods  are  still  effective  for  only  small  problems.  Reasonable 
computation  times  are  being  reported  for  networks  with  no  more  than  100 
activities  (  Talbot  and  Patterson  (1976).  Patterson  (1964)  ). 

9.42.  Heuristic  Methods 

Since  the  exact  methods  are  ineffective  for  larger  size  projects,  various 
heuristic  procedures  have  been  developed.  Heuristic  procedures  can  not 
guarantee  optimality,  but  they  are  able  to  handle  very  large  networks,  and.  in 
general,  may  give  good  enough  schedules  to  serve  the  operational  needs  of 
users. 

Heuristic  methods  are  schemes  for  assigning  activity  priorities  in  making 
the  activity  sequencing  decisions  required  for  resolution  of  resource  conflicts. 
There  are  a  large  number  of  heuristic  procedures,  some  of  which  can  be  found 
in  the  open  literature,  while  others  are  in  the  form  of  proprietary  commercial 
computer  software.  The  main  distinguishing  feature  among  heuristic  methods  is 
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the  measure  that  the  heuristic  uses  to  assign  the  priorities.  A  survey  by  Davis 
(1973)  is  given  for  the  heuristic  methods  in  the  open  literature. 

Heuristics  can  be  classified  into  parallel  or  serial  routines,  depending  upon 
whether  the  priorities  are  assigned  during  scheduling  or  before  any  activities 
are  scheduled,  respectively.  Although,  there  are  some  evidence  that  methods 
with  parallel  routines  perform  better  than  those  using  serial  routines,  there  is 
no  formal  study  that  proves  that  superior  methods  always  fall  into  one  of  the 
two  categories. 

A  list  of  heuristic  sequencing  rules  is  given  in  Table  3.1.  The  table  indicates 
the  researchers  who  found  the  rules  effective. 

A  study  by  Davis  and  Patterson  (1075)  provides  a  comparison  of  the  rela¬ 
tive  effectiveness  of  various  heuristic  methods.  The  study  covered  the  first 
seven  rules  in  Table  3.1.  Reported  results  show  that  the  Minimum  Slack  rule 
(also  called  the  Least  Total  Float  rule)  produced  the  best  results  for  all  perfor¬ 
mance  parameters.  The  Longest  Remaining  Path  and  Minimum  Late  Finish  rules 
performed  close  to  the  Minimum  Slack  rule.  In  a  more  recent  study,  Thesen 
(1976)  introduced  urgency  factors  to  activities,  which  are  calculated  as  linear 
combinations  of  activity  resource  utilization  rates  and  the  impacts  of  delays  in 
the  activity  start  times  on  the  schedule.  The  urgency  factors  are  used  to  find 
the  set  of  starting  activities  for  the  scheduling  period.  A  multi-dimensional 
knapsack  problem  which  maximizes  the  sum  of  urgency  factors  of  activities 
that  can  start  in  the  scheduling  period  within  resource  limits,  is  solved  at  each 
scheduling  period  to  find  the  starting  activities.  Thesen  also  reported  that  his 


Priority  Rule 

Explanation 

Researchers 

Minimum  Slack 

Priority  is  given  to  the  minimum  slack 
activity  A/in  [1$  -  £5*] 

Mize.  Fendley. 

Patterson.  Yiest 

Longest 

Remaining 

Path 

Priority  is  given  to  the  activity  with 
longest  remaining  series  of  activities. 
Calculations  are  made  using  pre¬ 
cedence  relations. 

Verhines,  Brand. 

Meyer,  and  Shaffer 

Minimum  Late 
Finish  Time 

Priority  is  given  to  the  activity  has 
minimum  late  finish  time,  ie.  as  calcu¬ 
lated  in  CPM  method. 

Pascoe.  Mueller- 

Mehrback.  and 

Gounguet 

Greatest 

Resource 

Demand 

Priority  is  given  to  the  activity  that  has 
maximum  total  resource 

requirements.ie.  =  flf<  £  a*, 

* 

Knight  and  Patterson 

Greatest 
Resource  Utili- 
zation 

Priority  is  given  to  the  activities  which 
results  in  maximum  resource  utilization 
at  each  scheduling  period.  This 
requires  0-1  integer  programming. 

RAMPS 

Shortest 

Imminent 

Operation 

Priority  is  given  to  the  activity  with 
shortest  duration 

Patterson 

Most  Jobs  Pos- 
sible 

Priority  is  given  to  the  combination  of 
activities  which  results  the  greatest 
number  of  activities  being  scheduled  in 
each  scheduling  period.  This  requires 
0-1  integer  programming. 

Some  commercial 

programs 

Hybrid  Rules 

Priority  is  given  to  the  activity  which 
has  highest  score  from  the  combina¬ 
tions  of  some  of  the  previous  rules 

Thesen 

TABLE  3.1.  Summary  of  Common  Priority  Rules  that  are  found  in  open  litera 
ture. 
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approach  was  superior  to  the  Minimum  Slack  priority  rule  according  to  his  com¬ 
parative  experiments. 

9.5.  Generalized  Resource  Constrained  Project  Scheduling  Problem 

In  this  section,  we  will  cover  the  methods  that  are  used  in  solving  the 
resource  constrained  problem  with  generalizations  of  some  aspects  of  the  com¬ 
mon  problem.  Sue  to  the  our  particular  interest  in  the  work  characteristics  in 
the  shipyard,  two  kinds  of  generalized  models  will  be  reviewed.  First,  the 
methods  that  relax  the  assumption  of  fixed  resource  application  rates  and  fixed 
durations  will  be  reviewed.  Second,  the  methods  that  handle  the  stochastic  ele¬ 
ments  in  project  scheduling  will  be  covered. 

3.5.1.  Hie  Methods  for  Resource  Application  in  the  Generalized  Problem 

The  first  two  published  models  which  allow  activities  alternative  rates  of 
resource  applications  were  Wiest's  SPAR1(1967)  model  (Scheduling  Program  for 
Allocation  of  Resources)  and  RAMPS  (Resource  Allocation  and  Multi  Project 
Scheduling)  by  Moshman  et  al.  (1963).  Both  procedures  consider  that  an 
activity  may  be  performed  either  at  a  normal  pace,  at  the  fastest  pace,  or  the 
slowest  pace  by  using  a  normai,  maximum,  or  minimum  amount  of  resources, 
respectively.  Although  the  RAMPS  software  has  been  widely  distributed,  the 
computational  steps  involved  are  nowhere  fully  described  in  detail. 

SPARl  basically  assumes  three  discrete  levels  of  resource  applications  for 
each  activity  (normal  crew  size,  maximum  crew  size,  and  minimum  crew  size). 
An  intermediate  level  can  be  applied  to  an  activity  only  in  special  cir- 
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cumstances.  The  detailed  procedure  of  the  SPAR1  model  will  be  covered  in 
Chapter  7. 

Weglarz  (1981)  studied  properties  of  optimum  solutions  for  doubly  con¬ 
strained  scheduling  problems  when  "performing  speeds”  of  activities  are  con¬ 
tinuous  functions  of  resource  application  In  his  model,  resources  are  not  only 
constrained  by  a  maximum  usage  rate  at  every  moment  but  also  by  a  maximum 
cumulative  resource  expenditure.  The  performing  speed-resource  usage  func¬ 
tion  for  each  activity  is  a  production  function  for  the  activity,  relating  the  rate 
of  completion  to  the  rate  of  resource  application.  Although  he  generalizes  the 
results  for  the  cases  when  activities  use  more  than  one  resource  or  have  upper 
and  lower  bounds  on  their  resource  usage,  it  is  not  clear  that  the  model  works 
in  the  multi-resource  case.  He  suggests  a  non-linear  programming  model  to 
find  the  optimal  solution  for  various  generalized  cases  of  the  resource  con¬ 
strained  problem  (  Veglarz  (1978),(1979)  ).  His  solution  technique  is  only 
applied  to  very  small  projects,  as  it  involves  modeling  and  computational 
difficulties  for  larger  size  problems. 

Talbot  (1982)  introduced  an  integer  programming  formulation  and  a  solu¬ 
tion  technique  for  the  case  that  each  activity  has  several  alternative  modes  of 
accomplishment  in  the  project.  Each  mode  has  a  different  magnitude  and  mix 
of  resource  usage.  Detailed  explanation  of  his  method  will  be  given  in  Chapter  7. 

Leachman  (1983)  considered  continuous  ranges  for  activity  intensity,  ie, 
an  index  of  the  rate  of  resource  applications  to  the  activity,  in  the  resource  lev¬ 
eling  problem.  This  continuous  intensity  range  is  adopted  in  the  approach 
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taken  in  Chapter  4  and  5. 

3.52.  The  Methods  that  Handle  the  Stochastic  Elements  In  the  Project  Schedul¬ 
ing  Problem 

Probabilistic  project  networking  represents  an  extension  of  classical  PERT 
in  order  to  incorporate  cases  in  which  the  way  the  project  continues  from  some 
point  on  may  depend  on  the  uncertain  outcome  of  one  or  several  activities. 
Such  situations  are  typical  with  R&D  projects.  One  of  the  best  known  systems  is 
GERT  (  Graphical  Evaluation  and  Review  Technique  )  originally  developed  by 
Pritsker  and  Happ  (1966).  Pritsker  and  Whitehouse  (1966).  and  Vhitehouse 
(1973). 

GERT  provides  several  flexibility  features.  It  allows  probabilistic  branching, 
ie,  the  need  for  execution  of  some  activities  may  be  uncertain  in  the  project 
network.  It  allows  looping  in  the  network,  ie,  depending  on  the  outcome  of  an 
uncertain  event,  a  branch  of  activities  that  has  been  executed  before  may  be 
undertaken  again.  It  allows  alternative  distributions  for  the  probabilistic 
activity  durations.  (Recall  that  the  PERT  method  assumes  beta  distribution  for 
all  the  activities  in  the  project.) 

Usually  GERT  network  analysis  is  carried  out  by  simulation  during  which 
statistics  on  the  project  times  and  cost  are  collected.  The  GERT  method  basi¬ 
cally  concentrates  on  time-cost  tradeoff  studies,  as  it  does  not  provide  any 
treatment  of  resource  constraints,  resource  application  flexibilities,  or  uncer¬ 


tain  resource  requirements  by  activities. 


VERT  (  Venture  Evaluation  and  Review  Technique  )  is  an  extension  of  GERT. 
It  provides  the  capability  to  model  decisions  within  the  network  in  terms  of 
times,  costs,  and  performances.  GERT  can  analyze  times  and  cost  only;  hence 
VERT  is  a  powerful  project  evaluation  technique  (  Mueller  (1972),  Mueller  and 


Digman  (1981)  ). 
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4.  MODELING  ASSUMPTIONS  OF  A  NEW  SCHEDULING  SYSTEM  FOR  THE  PROJECT- 
ORIENTED  PRODUCTION  ENVIRONMENT 

The  review  of  project  scheduling  methods  shows  that  no  one  of  the 
methods  alone  proposed  in  the  literature  meets  the  needs  arising  from  the 
characteristics  of  work  in  project-oriented  production  systems.  So,  a  new 
approach  of  planning  and  scheduling  the  project  in  such  systems  is  needed.  In 
this  chapter,  the  modeling  assumptions  for  the  new  system  will  be  put  forward. 

4.1.  A  Hierarchical  Approach  for  Multi-Project  Planning  and  Scheduling 

A  multi-project  environment  is  assumed  for  the  project  oriented  produc¬ 
tion  system.  In  a  multi-project  environment  where  multiple  projects  might  be 
executed  at  the  same  time,  one  could  construct  a  large  network  by  combining 
all  of  the  project  networks.  Scheduling  of  activities  would  be  done  from  the 
large  network,  utilizing  total  resources  available  in  the  production  system.  How- 
eve.-,  there  are  some  disadvantages  of  this  approach.  The  number  of  activities 
in  the  large  network  may  become  so  large  that  the  scheduling  problem  becomes 
intracable.  In  application,  each  project  has  different  management.  Scheduling 
of  activities  and  allocation  of  resources  between  activities  within  a  project  is 
decided  by  the  project's  management.  The  approach  introduced  above  requires 
scheduling  and  resource  allocation  decisions  to  be  made  by  only  one  manage¬ 
ment  for  all  projects.  This  may  create  some  managerial  difficulties. 

An  alternative  to  the  above  approach,  namely  a  hierarchical  project  plan¬ 
ning  approach,  can  be  used  (R.  Leachman  and  J.  Boysen  (1984).).  In  this 
approach,  first,  allocation  of  resources  is  made  between  projects.  Then,  pre- 


allocated  resources  to  a  project  are  used  as  capacities  of  resources  for  the  pro¬ 
ject  which  will  be  allocated  among  project  activities.  The  principal  advantages 
of  the  hierarchical  approach  is  its  suitability  to  organizational  structure  and 
that  it  makes  both  multi-project  planning  and  scheduling  of  individual  projects 
easier  to  handle.  Further  advantages  of  hierarchical  project  planning  and 
scheduling  approach  are  given  in  R.C.  Leach  man  and  J.  Boysen  (1964). 

4.2.  Project  Activities 

An  activity  is  a  part  of  the  project  which  has  a  beginning  point  and  an  end¬ 
ing  point  which  are  clearly  identifiable.  It  has  a  duration.  It  requires  resources 
of  different  types  for  its  execution. 

4.2.1.  Precedence  Relations  between  Activities 

We  shall  consider  the  case  when  precedence  relations  between  activities 
may  be  either  strict  or  overlapping.  Overlapping  activity  precedence  may  be 
used  in  reducing  the  network  size.  The  precedence  of  this  type  is  frequently 
useful  in  the  networks  modeling  ship  overhauls.  So.  project  networks  con¬ 
sidered  here  include  both  types  of  precedence  relations. 

4.2.2.  Flexible  Resource  Applications  for  Activities 

We  shall  assume  that  all  the  different  types  of  resources  which  are  required 
by  an  activity  are  applied  proportionally  through  the  activity  execution.  In 
other  words,  each  activity  utilizes  a  constant  mix  of  resources  as  it  progresses. 
When  the  rates  of  application  of  different  types  of  resources  to  the  activity  are 
proportional,  the  rates  can  be  indexed  in  terms  of  one  rate  which  will  be  called 
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the  intensity  of  the  activity.  The  rate  of  resource  application  of  an  activity 
(activity  intensity)  is  assumed  to  be  variable  within  upper  and  lower  limits. 
Each  activity  is  assumed  to  have  a  normal  resource  application  rate  and  a  nor¬ 
mal  duration  corresponding  to  constant  normal  intensity.  For  a  lower  rate  of 
resource  application,  activity  duration  will  be  longer  or,  for  a  higher  level  of 
resource  application,  activity  duration  will  be  shorter. 

The  flexibility  in  resource  application  to  activities  has  significant  impor¬ 
tance  in  project  scheduling.  The  scheduling  problem  in  the  case  of  fixed 
resource  application  rates  and  fixed  durations  is  just  to  decide  starting  times  of 
project  activities  which  will  not  violate  resource  constraints.  For  flexible 
resource  application  rates  to  activities,  the  problem  requires  deciding  resource 
application  rates  of  activities  in  addition  to  their  start  times. 

4.2.3.  Stochastic  Resource  Requirements  For  Activities 

Resource  requirements  of  an  activity  are  best  viewed  as  random  variables 
with  some  distributions  in  project-oriented  production  systems.  This  is  also  true 
for  construction  projects  but  the  projects  undertaken  by  construction  firms 
are  often  not  similar  to  one  another.  When  projects  are  different,  collection  of 
historical  resource  requirement  data  for  activities  is  not  possible.  Projects  are 
similar  or  even  identical  in  project-oriented  production  systems.  Projects  may 
be  executed  a  large  number  of  times  in  such  systems.  Consequently,  it  may  be 
possible  to  collect  data  necessary  for  finding  distributions  of  resource  require¬ 


ments  of  activities. 
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4.3.  Resources  Used  In  Projects 

In  general,  two  types  of  resources  are  used  in  projects.  These  are, 

a)  Non-storable  Resources 

Example:  Trade  labor  of  different  types,  machine  or  equipment  services 
used  in  production. 

b)  Storable  Resources 
Example:  Roney,  material,  etc. 

Our  main  concern  will  be  given  to  scarce  non-storable  resources  such  as 
trade  labor  from  different  shops.  Project-oriented  production  systems  may 
include  a  large  number  of  trade  labor  "shops'*.  Each  shop  consists  of  one  type 
of  skilled  trade  labor.  Employees  in  these  shops  work  on  a  fixed  salary  basis. 
The  cost  of  labor  to  the  production  system  is  mostly  fixed,  perhaps  with  the 
possibility  of  incremental  overtime. 

Each  shop  has  a  capacity  at  a  certain  time,  ie.  a  number  of  man-hours 
available  per  day.  Capacity  may  change  with  time  but  capacity  available  in  a 
given  time  period  is  inflexible.  The  total  shop  capacity  is  shared  by  the  ongoing 
projects  in  the  system.  The  aggregate  planning  level  allocates  the  shop  capaci¬ 
ties  between  projects  as  discussed  earlier.  Thus  each  project  has  preallocated 
capacity  from  all  shops.  This  capacity  is  to  be  allocated  between  activities  of 
the  project  that  are  ready  for  execution  at  that  time. 

Project  cost  for  the  production  system  largely  depends  on  project  time. 
The  objective  of  scheduling  in  such  a  system  is  to  finish  the  project  as  early  as 
possible  within  given  resource  capacities.  Equivalently,  the  objective  is  to 
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maximize  the  throughput  ,  ie,  the  number  of  projects  completed  by  the  system 
within  a  given  time  horizon. 

Storable  resources  may  be  limited  for  usage  per  period  as  non-storable 
resources  are,  and  also,  they  may  be  limited  in  their  cumulative  consumption 
over  the  period  of  project  duration.  The  latter  type  of  resources  is  of  secondary 
concern  in  our  case,  so  that  the  method  of  scheduling  is  first  designed  assum¬ 
ing  that  we  have  no  limitations  on  the  cumulative  consumption  of  resources. 
The  extension  of  the  method  with  necessary  adjustments  to  handle  cumulative 
resource  consumption  limitations  will  be  given  in  Appendix  2. 

4.4.  Stochastic  Project  Networks 

In  projects  like  overhauling  a  ship,  building  a  ship  or  an  airplane,  advanced 
testing  of  subsystems  is  needed  following  their  production  or  repair.  For  sub¬ 
systems  involving  high  safety  requirements,  eg.  nuclear  power  generating  sys¬ 
tem  in  nuclear  submarines,  testing  and  repair  rework  must  be  undertaken  until 
performance  of  the  subsystem  achieves  very  tight  tolerance  limits.  This  may 
take  a  number  of  test-rework  cycles.  The  number  of  test-rework  cycles  for  the 
same  test  in  different  examples  of  similar  projects  may  vary  greatly. 

We  refer  to  activities  involving  testing  of  subsystems  as  test  activities.  The 
outcome  of  a  test  activity  can  be  a  success  or  a  failure.  If  it  is  a  success  then  its 
regular  follower  activities  may  start.  If  it  is  a  failure,  depending  on  the  rework 
type  needed,  one  of  the  possible  small  subnetworks  of  rework  activities  must  be 
executed  before  the  test  activity  may  be  undertaken  again. 


Assuming  that  all  the  test  activities  will  be  successes  on  the  first  trials,  we 
have  a  basic  deterministic  network  structure.  After  each  test  there  are  a 
number  of  small  subnetworks,  one  of  which  is  to  be  executed  after  a  failure 
depending  on  the  type  of  the  failure.  Consequently,  a  special  stochastic  net* 
work  structure  is  assumed  for  the  projects. 
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5.  MODEL  WITH  VARIABLE  INTENSITIES  FOR  ACTIVITIES 

As  discussed  before,  there  is  considerable  flexibility  available  to  manage¬ 
ment  in  loading  of  manpower  resources  on  each  activity  of  a  project.  Most  pro¬ 
duction  activities  may  be  progressed  with  varying  intensity,  ie,  varying  levels  of 
manpower  application.  In  this  chapter,  we  develop  the  model  which  allows  vari¬ 
able  intensities  for  activities. 


5.1.  linear  Activity  Analysis  Model 

We  suppose  a  project  is  characterized  by  an  activity-on-node  critical  path 

analysis  network  of  N  activities,  denoted  by  A\  .  Ag .  Aa . A*.  The  project 

requires  utilization  of  K  infinitely  divisible,  non-storable  resources,  denoted  by 

Rx .  Rg .  Rs . Rk  •  Following  Shephard  and  Mak  (1982)  and  Leacbman 

(1963).  the  rates  of  application  of  various  resources  to  a  given  activity  A,  are 
assumed  to  be  proportional  and  constant  on  discrete  time  intervals 

(0,1],(1.2] . (t-l,t] . For  convenience,  we  refer  to  the  time  interval  (t-l,t]  as 

period  t.  The  rates  of  resource  applications  to  activity  A  in  period  t  are 
indexed  by  an  intensify  function  c((t),  taken  with  technical  coefficients  as  fol¬ 
lows: 


the  total  amount  of  exogenous  resource  k  that  is  required  by  A»  from  its 
start  to  the  end  ,  k-1,2 . K  . 

To  illustrate  this  model,  if  an  activity  A*  is  operated  continuously  at  a  con¬ 
stant  intensity  *t  .  then  its  duration  is  1/s*  .  and  resource  Rg  is  loaded  at  the 
rate  of  a*(*t  per  unit  time  over  the  duration.  Note  that  most  treatments  of 
resource  loading  based  on  ordinary  critical  path  analysis  implicitly  assume 
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activities  are  operated  in  this  constant  loading  fashion. 

We  shall  assume  each  activity  At  has  a  normal  duration  d*  which 
corresponds  to  normal  application  rates  of  resources  to  the  activity,  applied 
continuously  until  it  is  complete.  Then,  normal  intensity  is  defined  as, 

*.  =  £  <5i> 

With  this  assignment,  the  intensity  function  s<(i)  corresponds  to  the  frac¬ 
tion  of  the  activity  done  in  period  t.  So,  summation  of  st(f)  between  start  and 
finish  must  add  up  to  unity,  ie, 

t  **(0=1-  (5.2) 

!»«,♦! 

where  t,  *  tj  are  the  start  and  finish  times  of  A*. 

The  intensity  of  each  activity  is  assumed  to  have  upper  and  lower 
bounds  corresponding  to  practical  minimum  and  maximum  application  rates  of 
resources.  That  is.  from  the  start  of  the  activity  until  its  completion,  activity 
intensity  lies  in  the  range 

(5.3) 

where  the  bounds  jfc  and  It  are  givenf.  Before  the  start  and  after  the  comple¬ 
tion  of  Ai  ,  *i(t )  is  zero. 

5.2.  Intermediate  Activity  Transfers 

In  critical  path  analysis,  activities  are  usually  related  by  strict  precedence. 

Occasionally,  in  order  to  reduce  the  number  of  the  activities  required,  overlap 

t  In  the  case  that  resource  capacities  decrease  with  time,  activities  map  be  disrupted,  in 
which  case  activity  intensity  is  forced  to  aero  for  a  period  between  start  and  completion  of 
the  activity.  Disruption  of  activities  occurs  only  as  last  recourse  in  the  heuristic  algorithms 
to  be  developed. 
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relationships  may  be  specified.  (See  Elmaghraby  (1977).).  As  an  example,  sup¬ 
pose  A,  may  not  start  uniil  after  At  has  performed  a  fraction  of  its  work.  If 
<5*  is  the  start  time  of  A*,  and  Sj  is  the  start  time  of  Aj.  then  Sj  is  constrained  by 


^  +  Sii  s  (5.4) 

Note  that  /y  =  1  corresponds  to  strict  precedence  and  fqdi  is  approximated  to 

closest  integer  value. 

When  activity  intensities  are  variable,  the  overlap  relationship  (5.4)  must 
be  generalized.  To  do  this,  we  shall  use  the  concept  from  production  theory  of 
intermediate  product  transfer  between  dependent  (ie.  arc-connected  )  activi¬ 
ties.  as  set  forth  in  Shephard.  Al-Ayat  and  Leachman  (1977).  and  Leachman 
(1982). 


We  shall  say  that  there  is  a  progress  lag  of  between  A*  and  A,,  ie.  A» 
must  be  100 (/^+i)%  before  Ay  can  be  1005%  complete.  Formally,  given  an  /v 
satisfying  0</^sl  ,  we  have  ,  for  all  t. 


*t(T) 

if 


-/<,)  ^  0* £  *<(T)  <1*0 

J  rml 

£«t(r)=1.0 


(55) 


Once  again,  /y  =  1  corresponds  to  strict  precedence. 

In  the  case  that  intensity  must  be  held  constant  at  unity  from  start  to 
finish  of  each  activity.  (5.5)  holds  for  all  t  if  and  only  if  (5.4)  bolds  and  also. 

Si  +  d»  s  Sj  +  (1  -fa)  dj  (5.6) 

Note  that  if  ,  then  (5.4)  implies  (5.6).  However,  for  the  progress  lag 

here,  if  di>dj  ,  it  is  not  sufficient  to  simply  separate  the  start  times  of  constant 


intensity  activities  A%  and  Aj  according  to  the  lag;  the  activities  must  also  have 
finish  times  separated  according  to  the  lag.  as  expressed  in  (5.6).  When  inten¬ 
sity  is  allowed  to  vary,  the  more  general  expression  (5.5)  applies. 

When  Ai  and  Aj  are  related  by  a  progress  lag.  we  shall  say  there  is  a  flow 
transfer  between  them.  If  we  view  as  producing  intermediate  product 
required  as  input  by  Aj  then  after  an  initial  lag.  A<  continuously  supplies  inter¬ 
mediate  product  to  Aj,  allowing  progress  to  be  made  towards  completion  of  Aj. 
When  fij-l,  there  is  a  strict  precedence  between  A^  and  Aj  ,ie  .the  transfer  of 
intermediate  product  is  entirely  event-based.  We  shall  assume  in  the  following 
that  the  activity  network  consists  of  both  event-based  and  flow  transfers  of 


intermediate  products. 
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0.  ALGORITHMS  TOR  INTENSITY  ASSIGNMENT 

6.1.  Introduction 

This  chapter  presents  two  algorithms  (Upgrading- Only  and  Upgrading  and 
Downgrading)  for  resource  constrained  project  scheduling  with  variable  intensi¬ 
ties  for  activities.  An  intensity  assignment  algorithm  can  be  visualized  as  a 
decision  maker  which  decides  the  intensities  of  activities  which  will  be  main¬ 
tained  until  the  next  decision  time.  Through  intensity  assignment  to  an  activity, 
resourcesf  used  by  this  activity  are  also  allocated,  as  explained  in  section  5.1. 
At  each  decision  time,  the  algorithms  allocate  the  resources  up  to  their  capaci¬ 
ties  among  activities  that  are  already  started  or  ready  to  start  at  that  time. 
Decisions  have  to  be  revised  at  time  points  when  there  is  a  change  in  the  system 
that  may  effect  the  intensity  assignment.  These  changes  are  completion  of  an 
activity,  the  event  of  an  activity  becoming  ready  to  start,  and  a  change  in  capa¬ 
cities  of  resources.  Tor  the  moment,  it  is  assumed  that  all  resource  capacities 
are  constant.  This  assumption  will  be  relaxed  in  section  5.5. 

The  main  objective  of  the  scheduling  problem  is  to  find  a  schedule  which 
gives  the  least  project  completion  time  within  the  capacity  restrictions  and 
event  due  dates  which  are  determined  by  the  aggregate  plan.  The  intensity 
assignment  algorithms  to  be  developed  consist  of  ’clever’  heuristic  decision 
rules  to  be  applied  at  each  decision  time.  Consistent  application  of  these 
heuristics  will  result  in  a  schedule  with  a  'good'  project  finish  time. 


t  Only  non-storable  resources  ere  considered  in  this  chapter.  The  eztenaion  of  the  metho-' 
dolcgy  introduced  in  this  chapter  for  storable  resources  will  be  given  in  Appendix  A2. 


The  basic  intuition  behind  the  logic  of  both  algorithms  is  equalizing  rela¬ 
tive  lateness  (or  earliness)  of  activities  against  their  target  finish  times  through 
assigning  appropriate  intensities  to  the  activities,  argued  as  follows. 

A  Late  Finish  Time  of  an  activity  can  be  visualized  as  a  due  date  for  the 
activity  finish,  if  it  is  not  specified  by  the  aggregate  plan,  such  that  if  the 
activity  finishes  any  later  than  the  due  date,  it  will  lengthen  the  project  com¬ 
pletion  time  or  violate  the  aggregate  plan.  Since  this  is  true  for  all  activities, 
the  decrease  in  project  finish  time  can  only  be  achieved  if  all  activities  finish 
early  with  respect  to  their  due  dates.  Any  one  activity  being  late  with  respect 
to  its  due  date  is  enough  to  increase  the  project  finish  even  if  all  the  others  are 
early  with  respect  to  their  due  dates.  The  sensible  objective  at  each  decision 
time  is  to  make  all  activities  ,  if  possible,  to  be  as  early  as  possible  with  respect 
to  their  target  late  finish  times.  This  is  accomplished  by  equalizing  earliness. 

6.1.1.  Activity  Target  Finish  Time 

The  target  finish  time  of  an  activity  is  taken  as  the  Late  Finish  Time  of  the 
activity  as  it  is  computed  using  normal  durations  in  ordinary  critical  path  cal¬ 
culations  subject  to  the  due  dates  specified  by  the  aggregate  plan.  Before  each 
intensity  assignment,  the  critical  path  Late  Finish  Time  calculation  is  per¬ 
formed  for  the  unexecuted  part  of  the  activity  network,  ie,  the  activities  that 
are  not  started  yet  and  the  activities  that  are  started  but  not  finished  yet.  In 
this  calculation,  it  is  assumed  that  activities  will  use  constant  normal  intensi¬ 
ties  between  their  start  and  finish  times,  ie.  the  same  assumption  as  in  critical 
path  analysis.  For  the  activities  that  are  already  started  and  already  have 
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intensities  assigned,  normal  intensity  is  assumed  for  their  uncompleted  frac¬ 
tions. 


6.12.  Notation 

We  shall  use  the  following  notation  (in  addition  to  notation  set  forth  in 

Chapter  5). 

A*  set  of  indices  of  activities  to  which  At  supplies  intermediate  product  ( 
followers  ). 

r<  set  of  indices  of  activities  supplying  At  with  intermediate  product 
(predecessors). 

t  current  time  period,  ie.  reallocation  of  resources  is  considered  at  time 

(*-!>• 

*i(f )  intensity  of  At  in  time  period  t  . 

<5«  actual  starting  time  of  A%. 

Ft  actual  finishing  time  period  of  At . 

Ft(t )  finishing  time  of  At  projected  at  time  t  . 

FTtjt{t)  flow  transfer  time  for  follower  At  tie,  the  time  that  activity  i  will  be  pro¬ 
gressed  enough  to  allow  its  immediate  follower  l  to  start,  projected  at 
time  t  . 

pcvt(t )  fraction  of  At  completed  at  time  t  .  (  s  £  *.(r)) 


LFt(t)  target  late  finish  time  of  At  computed  at  time  t  .  computed  by  an  ordi¬ 
nary  (ie,  resource-unconstrained  )  critical  path  analysis  backward  pass 


■v.  ---/v  w -v .*•: -vt- 
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assuming  normal  intensities  for  all  activities.  Detailed  calculation  of 
LFx(t)  will  be  described  in  section  6.1.3. 

Isci(t)  score  of  lateness  of  Ai  at  time  t  ,ie,  lateness  with  respect  to  its  target 
late  finish  time. 

esc»(t)  score  of  earliness  of  yi,  at  time  t  ,ie,  earliness  with  respect  to  its  target 
late  finish  time. 

rM(t )  capacity  of  resource  R*  in  time  period  t  . 

OTffit)  amount  of  resource  /4  available,  ie,  unallocated,  in  time  period  t  . 

IACTIVE  set  of  indices  of  currently  active  activities,  ie. 

|  0.  <1.0 

1READY  set  of  indices  of  activities  that  are  currently  ready  to  start  ie,  those 
activities  not  yet  active  whose  predecessors  are  progressed  enough  to 
satisfy  the  progress  lag  (equation  5.5). 

6.1.3.  Late  Finish  Updates 

Recalculations  of  late  finish  times  of  activities  are  made  before  every 
assignment  of  intensities,  that  is,  there  is  a  recalculation  of  the  times  used  as 
targets  for  activities  for  intensity  assignment  purposes.  Recall  that  parallel 
routines  update  the  priorities  in  the  same  way  for  the  fixed  intensity  case.  In 
calculating  these  times,  it  is  assumed  that  activities  which  are  not  finished  at 
the  time  of  update,  will  run  with  normal  intensity  from  that  time  on.  This  is  true 
for  active  activities,  regardless  of  current  intensity  level,  and  also  for  unstarted 


activities.  The  following  steps  are  performed  at  each  late  finish  update. 


1.  Find  Adjusted  Finish  Times  for  all  active  and  completed  activities  (Ai).  The 
projected  finish  time  Ft(t-1)  reflects  the  current  intensity  of  A*  ,ie. 
•i(f-l).  AFi(t-l)  is  the  adjusted  finish  time  of  A*  assuming  that  it  will 
have  unit  intensity  after  time  (< — 1).  It  is  calculated  as  follows: 

(a)  If  an  activity  is  finished  at  or  before  the  current  time  t-1  .ie. 
f  — 1  St  Ft,  then 

AF^t-l)  =  F, 

(b)  If  the  activity  is  in  progress,  then 
AF^t-i)  =  f-1  +  (  1.0 -pcv^t-l)  )  di 

(c)  If  the  activity  is  not  started  yet,  then  set  the  adjusted  finish  time  to  a 
big  number  ,  T  ,  so  that  it  is  distinguished  in  future  calculations. 
AFi(t~l)  =  T  if  t-1  <  Si 

l.  Run  CPM  type  calculations  to  find  late  finish  times  as  follows: 

(a)  Set  the  early  finish  times  of  activities  that  are  finished  or  started  but 
not  finished  at  time  (t  —  1)- 

EFi(t-l)  -  AFiit-1)  if  AFi(t-l)*T 

where  EFiit  -1)  is  early  finish  time  of  At  at  time  t*l  . 

(b)  Make  forward  calculations  to  find  Early  Finish  times  for  unstarted 


activities,  ie,  for  and  A%  unstarted, 

EFi  -  Ufut  [EFj]  4  dt 


(c)  Make  backwards  pass  to  find  target  late  finish  time  (  LFt(t )  ). 


(i)  For  each  activity  A*  that  has  a  due  date  ( DD *)  specified  by  the 
aggregate  plan,  set  its  Late  finish  time  equal  to  the  due  date. 

(It  is  assumed  that  all  activities  with  no  followers  have  due  dates 
specified  by  the  aggregate  plan.  Other  activities  may  have  due 
dates  as  well.) 

(ii)  Make  backward  calculations  for  the  activities  that  have  fol¬ 
lowers.  For  such  that  A.  is  nonempty,  set 

LFAt)  =  Jfin  [LFj  -  dj) 


6.1.4.  Equalising  Relative  Lateness  (or  Earliness) 

As  discussed  before,  the  ideal  objective  of  intensity  assignment  is  equaliz¬ 
ing  earliness.  The  problem  can  be  formulated  mathematically  as  follows. 

Let  S=i4,  ,  •  •  •  ,  Ai  ,  •  •  •  ,Ai  be  the  set  of  activities  that  are  candidates 
for  intensity  assignments  at  time  t  -1  . 

Then  we  wish  to 

Uin  £  (Fi(t)-m(t- D)e 

Id 

subject  to: 


Fi<0  = 


_  (  1  -  pcTj(f -1)  ) 


«i(0 


♦  (1-1) 


for  all  Ai  t  S  (  Finish  Time  ) 


£  o*i  *i(0<  rk(t) 


tor  all  k  (Resource  Availablity) 


pcvt(t -l)+zl(t)*pcvk(t -l)+**(0-/«  for  kchi.  AitS,  A*tS  (Flow  transfers) 
li  <  et(t)  £  Ft  or  tt(t )  =  o  for  all  Ai  t  S  (Intensity  bounds) 

Note  :  Ci(f)‘s  and  Fi(t)' s  are  decision  variables.  Other  parameters  have 
known  values  at  time  t. 

The  above  formulation  has  a  quadratic  objective  function  with  some 
integer  variables  used  to  handle  the  complication  caused  by  the  fact  that  **(0 
has  to  be  either  within  its  upper  and  lower  bounds  or  zero.  The  problem  does 
not  have  a  structure  that  admits  an  easy  solution  methodology  to  find  the 
optimal  solution. 

A  method  for  a  simplified  version  of  this  problem  is  solved  by  Bitran  and 
Has  (1981)  for  the  problem  of  equalising  run-out  times  of  inventories  of  the 
items  in  the  same  product  family  in  a  hierarchical  production  planning  system 
In  their  model,  there  is  only  one  capacity  constraint  and  decision  variables 
have  only  upper  and  lower  bounds.  Their  method  is  not  applicable  to  our  prob¬ 
lem  with  multiple  capacity  constraints.  Consequently,  we  shall  develop  a  greedy 
priority  approach  to  allocate  resources  to  meet  target  finish  times  as  discussed 
below. 

0.2.  Greedy  Method  of  Assigning  Intensities  to  Activities. 

At  each  decision  time  a  list  of  activities  which  are  candidates  for  intensity 
assignment  is  generated.  The  relative  lateness  considering  target  Late  Finish 
Times  of  the  activities  on  the  list  is  assessed.  Activities  on  the  list  are  priori¬ 
tized  using  lateness  scores  equal  to  (Ft(t)-LFt(t)),  whereby  the  larger  the  late¬ 
ness  score,  the  higher  the  priority.  These  priorities  determine  the  order  that 


activities  will  be  considered  for  intensity  assignment.  Note  that,  the  priority 
rule  selected  is  similar  to  the  'Minimum  Slack'  and  *Earliest  Late  Finish*  rules 
tested  and  found  superior  (in  the  fixed  intensity  case)  to  other  rules  by  Davis 
and  Patterson  (1975),  and  othersf. 

Three  different  passes  are  made  before  finalizing  the  intensity  assignments 
to  the  activities  in  the  candidate  list.  In  the  first  pass,  using  the  priority  list, 
activities  are  considered  one  by  one  until  all  of  them  are  considered.  For  each 
activity,  an  attempt  is  made  to  upgrade  the  intensity  of  the  activity  (or  start 
the  activity),  if  necessary,  up  to  the  critical  intensity  level  that  will  allow  the 
activity  to  just  meet  its  target  finish  time.  Note  that  resource  availabilities,  the 
intensity  upper  bound,  or  progress  lag  limitations  allowed  by  flow  transfers  may 
be  more  restricting  than  the  critical  intensity  level.  If  this  is  the  case,  the  most 
restricting  level  is  taken.  If  the  intensity  level  taken  is  under  the  lower  Intensity 
bound  of  the  activity  than  zero  intensity  is  assigned  to  the  activity,  ie,  the 
activity  remains  unscheduled,  at  this  pass. 

In  the  second  pass,  a  look  ahead-procedure  is  used  to  correct  mistakes 
that  may  be  made  in  earlier  intensity  assignments.  The  procedure  is  similar  to 
the  one  used  in  Wiest's  SPARl  model.  In  this  pass,  attempts  are  made  to 
upgrade  intensities  of  the  critical  and  not  yet  started  activities  by  reschedul¬ 
ing  (ie,  delaying  )  the  truly  alack  ftt  activities  started  in  earlier  periods.  The 

only  reason  a  critical  activity  would  not  be  started  at  the  first  pass  is  if 
?  Iii2«  (1964),  Pascoc  (1965),  Go&ruet  (1969),  Patterson  (1975). 

tt  An  activity  is  called  critical  if  it  can  not  catch  up  to  its  late  finish  time  running  with  nor¬ 
mal  intensity  from  the  current  time  t.  ie,  its  lateness  score  is  negative, 
ttt  An  active  activity  is  called  truly  slack  if  by  starting  the  activity  in  period  (t+1)  with  nor¬ 
mal  intensity,  it  would  be  completed  before  its  target  late  finish  time. 
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resource  availability  for  the  activity  (or  its  predecessors)  is  restricting.  By 
rescheduling  the  truly  slack  activities,  the  resources  which  were  tentatively 
allocated  to  these  activities  are  made  available  to  such  critical  activities. 

in  the  third  pass,  using  the  priority  list,  activities  are  again  considered  one 
by  one  until  all  activities  are  considered.  For  each  activity,  an  attempt  is  made 
to  upgrade  the  intensity  of  the  activity  up  to  Its  upper  intensity  bound.  Again 
resource  availability  or  progress  lag  limitations  of  flow  transfers  may  be  res¬ 
tricting.  The  most  restrictive  level  of  intensity  is  assigned  to  the  activity,  if  this 
level  is  higher  than  the  intensity  lower  bound  of  the  activity.  Otherwise,  zero 
intensity  is  assigned  to  the  activity. 

Intensities  assigned  to  activities  at  the  third  pass  are  the  levels  that  activi¬ 
ties  will  operate  at  until  next  intensity  assignment  (ie,  decision)  time.  If 
resources  are  tight,  activities  which  are  at  the  top  of  the  priority  list  will  be 
allocated  resources,  while  activities  which  down  in  the  list  will  be  postponed. 
Eventually,  all  activities  so  postponed  become  critical  and  move  to  the  top  of 
the  priority  list. 

Using  this  basic  approach,  two  algorithms  which  reflect  two  different 
management  policies  concerning  resource  reallocation  are  introduced  in  sec¬ 


tions  8.3  and  8.4. 


49 


6.3.  Upgrading-Only  Algorithm 

We  first  consider  the  case  where  management  does  not  allow  reduction  in 
the  level  of  resources  which  is  already  committed  to  an  activity  once  it  is 
started,  although  it  allows  increases  in  this  level  while  the  activity  is  in  pro¬ 
gress.  The  Upgrading-Only  algorithm  is  designed  to  cope  with  this  case.  This 
algorithm  reassigns  the  resources  to  the  activities  by  reassigning  intensities  to 
the  activities,  and,  it  is  employed  whenever  there  is  an  event  that  may  cause  a 
change  in  current  resource  allocation.  The  algorithm  is  initialized  as  follows. 

INITIALIZATION: 

MOTIVE  =  * 

IREADY  *  Ji  ]  T,  =  #J  for  all  i 

zt{t)  =  0.  . pcv,(f )  =  0,  for  all  t. 
ar*(l)  =  r*(!)  tor  all  k- 

The  first  event  is  the  start  of  the  project,  time  f=0  .  Follow-on  events  are 
calculated  by  the  algorithm.  When  an  event  that  may  cause  resource  realloca¬ 
tion  occurs,  the  algorithm  performs  the  following  steps: 

STEP-1 

Update  target  late  finish  times  LFi(t  - 1)  by  using  the  procedure 
described  in  section  6.1.3. 

Note:  time  (t-1)  denotes  the  time  just  before  resource  allocations,  and 
period  t  denotes  the  period  (f  —  l,t )  in  which  allocations  are  applied 


STEP-2  Upgrade  if  necessary  intensities  of  active  and  reedy  activities  so  that 
they  may  catch  up  to  their  late  finish  times  (targets),  as  follows. 


(a)  Update  Lateness  Scores  assuming  ready  activities  will  run  at  normal 
intensity. 

lsct(t  -1)  =  Fdt  -1)  -  LFiit  -1)  for  all  <  t  IACTIVE 
lsct{t  -1)  =  (t  -1)  +  di  -  LFi(t  -1)  for  all  i  e  IREADY 


(b)  Prioritize  activities  in  combined  list  ( IREADY  and  IACTIVE)  by  sorting 
them  in  decreasing  order  with  respect  to  Lateness  Scores.  In  the  case  of 
ties,  any  order  suffices  which  is  consistent  with  network  logic. 

(c)  Upgrade  intensities  of  activities  in  combined  list  one  by  one  according  to 
priority  as  follows: 

set. 


Vi(t)  =  Uin 


Intensity  Upper  Bound 


min  I  —  -— (f -1)  Resource  Availiblity 

*  l  J 

(1.0-pcvi(f-l)) 


1.0 


min 

i'Tt 


LF^t-U-t 

(1.0 -*>«,<(* -!)-/><) 


Target  Finish  Time 

1.0 


Flow  Transfer T 


Note  that  y«(f)  is  a  temporary  variable  to  calculate  s*(f)  .which  is  defined  in 
terms  of  yt(f )  as  follows: 


*,«)  = 


Vi(t)  if 
*<(f- 1)  if 


0  if 


yt(f)fcet(f-l)  and  y* (0*1* 
yi{t)<si{t- 1)  and  *i(f-l)fcfc 
V*(f  )<*x(t-l)  and 


t  Derivation  of  this  constraint  is  given  in  Appendix  A1 


Then  set 


Fiit)  =  t- 1  +  (1.0-pcv<(f  -1)) 


and  reset 

*rk{t)  «-  ark (t )  -  a *<(*<(1  )-*<(l-l)) 

Repeat  the  above  for  <  e  IACTIVE  U  READY  with  next  highest  priority,  until  all 
the  activities  in  combined  list  are  processed.  Then,  continue  to  step  3  . 

STEP-3  Try  to  upgrade  intensities  of  the  critical  and  not  yet  started  activities 
by  rescheduling  the  truly  slack  activities  started  in  earlier  periods.  The 
only  reason  a  critical  activity  would  not  be  started  at  step- 2  is  if  the 
available  resources  are  not  enough  to  start  the  activity  with  an  inten¬ 
sity  level  which  is  higher  than  its  lower  intensity  bound.  By  reschedul¬ 
ing  the  truly  slack  activities,  the  resources  which  are  used  by  these 
activities  are  made  available  to  such  critical  activities. 

This  is  done  as  follows: 

(a)  Find  a  critical  but  not  started  activity  from  ready  activity  list  (READY).  If 
there  is  not  such  an  activity,  or  if  we  have  attempted  all  to  upgrade  all 
such  activities  continue  to  step-4.  Otherwise,  let  At  be  such  an  activity. 
Then  it  satisfies 

isc<(l)  £  0  and  *<(1-1)  =  0  where  <  e  READY 


(b)  Determine  if  enough  resources  will  be  obtained  by  rescheduling  truly 
slack  activities  that  are  already  started. 
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(i)  Initialize  the  search  variables: 

meed*  =  a**  &  -  ar*(f)  for  all  k 
JRESCHED=  i 

Note  that  meed*  and  IRESCHED  are  temporary  variables,  meed* 
represents  further  resource  needed  from  resource  type  k  to  start  4  . 
IRESCHED  represents  the  list  of  activities  that  are  to  be  rescheduled  in 
order  to  start  A t . 

(ii)  Find  a  truly  slack  activity  from  the  active  activity  list  {IACTD'E).  Let  Aj 
be  such  an  activity.  Then  it  satisfies: 

f  +  di  sc  LFj 

If  there  is  such  an  Aj  .  continue  to  (iii). 

Otherwise,  go  to  (a)  to  get  another  critical  activity. 

(iii)  Check  if  the  resources  needed  by  At  have  been  obtained. 

Update: 

meed*  «-  meed*  -  akJ *j(t -1)  for  all  k 
then. 

If  meed*  <  0  for  all  k  then  continue  to  (c) 

Otherwise,  put  Aj  into  the  reschedule  list. 

IRESCHED  -  IRESCHED  +  ; 

Go  to  (ii)  to  find  another  activity  to  reschedule. 
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(c)  Reschedule  the  activities  in  the  list  1RESCHED  since  they  provide  enough 
resources  to  start  Ax.  and  upgrade  the  intensity  of  Ai- 

(i)  Update: 

ark(t )  *-  ark(t )  +  akj  *j(t  -1)  for  all  j  e  IRESCHED 
Mj(t)  =  0  for  all  j  t  IRESCHED 

pcvj  (0  =  0  for  all  j  t  IRESCHED 

1ACTIVE  •-  IACTIVE  -  j 
IREADY  -  IREADY  +  j 

(ii)  Upgrade  the  intensity  of  >1*.  This  can  be  done  by  following  the  same 
procedure  as  in  step*?  (c). 

(iii)  Go  to  (a)  to  find  another  critical  activity. 

STEP-4  Use  remaining  resources  for  further  upgrading  of  activities  in  com¬ 
bined  list. 


(a)  repeat  step  1(a) 

(b)  repeat  step  1(b) 


(c)  Upgrade  Intensities  using  same  rule  in  1(c)  except  revise  y<(0  as  follows: 


Vt(0  = 


min 

* 


ark(t) 


0*1 


Intensity  Upper  Bound 
+  **(£-1)  Resource  AvaUiblity 


nun  (l.O-pcViit-l)-/^)  — ^ now  7Vans/er 
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STEP-5  Make  necessary  updates: 

(a)  Record  the  start  of  each  new  activity:  if  A+  is  such  that  ej(<-l)  =  0  and 
**(/ )fcO  then, 

IREADY  =  IREADY  -  i 
I  ACTIVE  -  IACTJVE  +  i 

(b)  Update  flow  transfer  times  of  upgraded  activities 

FTu(t)  =  (/«  - pcvt (f-1))  for  it  r* 

(c)  Find  next  resource  allocation  time  fn: 


tn  =  Min 


FiV)  .  FTu(t) 


for  itIACTIVE  and  all  It 


Note:  resource  allocation  is  needed  either  when  an  activity  finishes  or 
when  new  activities  become  ready  to  start.  If  calculated  tn  is  decimal,  it 
will  be  rounded.  As  discussed  in  section  6.5.2,  if  there  is  a  capacity  change 
before  time  tn,  the  time  of  capacity  change  is  treated  as  the  next  event 

• 

(d)  Update  intensities  and  cumulative  intensities  during  (t,tn) 

*<(*)  =  *<(0  •  t  —t  +  i,...,fn 

pCVt(T)  =pCVi( T-l)  +  *i(t)  .  T=f....tn 

Fi(r)  =  /;(/).  T=f  +  1 . tn 

/T(u)(t)  =  FTu(t).  r=f  +  l . fB 

(e)  Update  ready  activity  list: 


1 

I 


I  READY  *-  IREADY+ j  for  all  j  satisfying  pcvj(fn)fc/y 

for  all  l  cKj 

(f)  Update  available  resources  to  reflect  activities  finishing  at  t 

07*£  (tn )  -  ark(t)  +  )  if  is  such  that 

pcvi(tn)=l.  and pcvt(t)>  1. 

(g)  Update  current  resource  allocation  timer 
t  -  tn 

STEP— 6  Check  if  all  activities  are  finished: 

If  IREADY - IA CTIVE-  ♦ ,  then  all  activities  are  done, 
ie.  end  of  scheduling. 

Otherwise  go  to  step  0  for  resource  reallocation. 
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6.4.  Upgrading  and  Downgrading  Algorithm. 

We  now  consider  the  case  where  management  allows  both  increase  and 
decrease  in  the  level  of  resources  which  is  assigned  to  an  activity  that  is  in  pro¬ 
gress.  The  Upgrading  and  Downgrading  algorithm  is  designed  to  cope  with  this 
case. 

The  Upgrading-Only  algorithm  is  easily  modified  for  this  case.  We  need  to 
perform  a  step  before  the  Upgrading-Only  algorithm  which  is  as  follows. 

STEP-0 

Set  (tentatively)  the  intensities  of  i  e  MOTIVE  to  their  lower  intensity 
bounds,  as  follows: 


«i(0=Ai  tor  it  MOTIVE 

=  «-l)  +  (1.0  -pmt(t-l))  ^2. 

Mi 

and  reset 

arg(t)-ark(t)  +  o«(*t(f  )-»*(< -1)) 

Then  perform  the  steps  of  the  Upgrading-Only  algorithm. 

As  can  be  seen,  all  active  activities  are  tentatively  set  to  lower  intensity 
bounds  in  order  to  increase  pools  of  available  resources  to  maximum  levels. 
The  activities  are  then  upgraded  according  to  urgency  in  steps  1  through  6  of 


the  Upgrading-Only  algorithm. 


6.5.  Handling  Time-Varying  Capacities  For  Resources. 

6.5.1.  Basic  approach 

In  project-oriented  production  systems  resource  capacities  available  to  a 
particular  project  change  frequently.  A  change  in  resource  capacities  is  one  of 
the  events  that  may  require  changes  in  the  intensities  of  the  activities  assigned 
using  one  of  the  algorithms.  If  at  a  decision  time  there  are  increases  in  the 
capacities  of  certain  types  of  resources,  the  intensity  assignment  algorithms 
would  be  used  to  reassign  the  intensities  of  the  activities  in  order  to  utilize  the 
extra  resources  available  from  the  capacity  change.  If  there  are  decreases  in 
capacities  of  some  resource  types,  withdrawal  of  resources  from  some  activities 
may  be  necessary.  Even  stopping  the  executions  of  some  activities  may  be 
necessary  but  stopping  activity  execution  is  considered  only  as  a  last  resort. 

A  heuristic  procedure  is  applied  in  each  decision  time  with  capacity 
changes.  The  heuristic  tries  to  restore  the  feasibility  from  infeasibilities  which 
may  be  caused  by  resource  capacity  changes.  If  there  is  no  infeasibility,  ie,  new 
capacities  for  the  resources  are  greater  than  the  total  resources  used  by  all 
activities  which  are  in  execution,  then  intensity  assignment  algorithms  are  used 
to  reallocate  the  resources  between  activities.  If  there  are  infeasibilities,  ie.  new 
capacities  of  some  of  the  resources  are  smaller  than  the  total  resources  used 
by  the  activities  which  are  currently  in  execution,  then  attempts  to  restore 
feasibility  will  be  made. 

A  priority  list  is  generated  for  the  activities  that  are  active  using  their  ear¬ 
liness  scores,  ie,  the  negatives  of  their  lateness  scores.  First,  using  this  priority 
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list,  activities  are  considered  for  intensity  decrease,  one  by  one,  until  all 
infeasibilities  are  eliminated  or  until  all  activities  are  considered.  Intensity  of 
each  activity  is  set  to  its  lower  bound  or  to  an  intensity  level  that  will  be  just 
enough  to  restore  feasibility.  Note  that  if  an  activity  does  not  use  the  resource 
types  that  are  causing  infeasibility  its  intensity  remains  same.  If  all  infeasibili¬ 
ties  are  not  eliminated  in  this  step,  then,  using  the  same  priority  list,  activities 
are  considered  for  disruption  until  feasibility  is  restored. 

6.5.2.  Capacity  Change  Procedure 

A  time  of  capacity  change  is  treated  as  an  event.  We  now  designate  capaci¬ 
ties  as  rk(t) .  k=l,...K  .  Suppose  for  some  k  .  rM(t-l)  *  rk(t )  ,ie  .  there  is  a  capa¬ 
city  change  event  at  time  t-1.  For  such  an  event,  the  following  steps  ore  per¬ 
formed  by  the  algorithm  before  the  steps  of  the  intensity  assignments  at  that 
event.  Let  Ar*(t)  be  amount  of  capacity  change  for  resource  type  k  at  time  t,  ie, 

STEP— 1  Update  following  parameter, 

or*  (* )  «-  arfc  (< )  +  Ar*(f )  all  k 

STEP— 2  Check  if  there  are  any  capacity  violations 

If  art(l)  2  0.  then  Go  to  5 
Otherwise  Go  to  3 

STEP-3  In  this  step,  we  downgrade  intensities  in  an  attempt  to  rectify  capacity 


violations. 


(a)  Update  target  late  finish  times  LFi(t-l)  by  using  the  procedure  described 
in  section  6.1.3. 

(b)  Calculate  Earliness  Scores  for  all  active  activities. 


asciit-l)  =  L/;(f-l)  fori  t  I  ACTIVE 

(c)  Prioritize  activities  by  sorting  active  activities  in  decreasing  order  with 
respect  to  Earliness  Scores.  In  the  case  of  a  tie  any  order  is  acceptable 
which  is  consistent  with  network  logic. 

(d)  Downgrade  intensities  of  active  activities  one  by  one  according  to  priority 
as  follows.  Choose  the  active  activity  with  the  highest  earliness  score,  and 
set 


*i(0  = 


Pi(t-l)  ifoki=0  for  off  k  such,  that  art(t)<0 

otherwise 


Max 

z*(f  -1  )-Aftn 

-art(f)' 

.  Zi 

* 

Update, 


Ft(t-l)  =  t-1  +  (1.0  -pei^f-l)) 
o r*(f)«-ar*(0  +  aM(zi(t  )—**(<  — l)) 


If  or*(f)iO  for  all  k  .then  ,  go  to  step  5  :  otherwise  repeat  the  above  rou¬ 
tine  for  each  itIACTTVE  in  order  of  highest  priority.  If  all  active  activities 
have  been  processed,  continue  to  step  4. 

STEP-4  In  this  step,  activities  are  disrupted  until  capacity  violations  are  elim¬ 
inated.  Retaining  the  prioritization  from  step  3,  activities  are  dis¬ 


rupted  one  by  one  according  to  priority,  as  follows. 


*t(f  -1)  if  aMsO  for  all  k  such  that  a r*(f )  >  0 
*»(f)  -  |  o.  otherwise 

Update, 

or*(f)  «-  arM(t)  +  a*  (*,(!)-*,  (l-i)) 

/<(f-l)  -  T  (T  is  a  big  number  ) 

If  **(0  *  0  .  go  to  step  5,  otherwise  repeat  the  above  for  itIACTIVE 
with  next  highest  priority. 

5  Go  to  Intensity  Assignment  Algorithm 


7.  VALIDATION  OF  INTENSITY  ASSIGNMENT  ALGORITHMS 

The  aim  of  this  chapter  is  to  provide  some  evidence  of  validity  of  the  inten¬ 
sity  assignment  algorithms  which  are  introduced  in  Chapter  6.  The  best  way  of 
doing  this  would  be  to  compare  the  schedules  that  are  obtained  from  intensity 
assignment  algorithms  to  the  optimal  schedule,  if  there  was  a  method  available 
which  provides  an  optimal  solution  to  the  problem.  Unfortunately,  such  a 
method  does  not  exist.  Another  way  would  be  to  compare  the  schedules  that 
are  obtained  from  intensity  assignment  algorithms  to  the  schedules  that  are 
obtained  from  some  other  methods  which  solve  the  problem  with  the  same 
assumptions  that  are  made  for  intensity  assignment  algorithms.  Unfortunately, 
no  methods,  that  we  know,  have  been  offered  in  the  literature  that  solve  the 
problem  with  exactly  the  same  assumptions.  But.  two  methods  which  solve  the 
problem  with  slightly  different  assumptions  are  available.  These  are  Wiest's 
heuristic  method  (SPARl)(l9B7),  and  Talbot's  implicit  enumeration  method 
(19B2).  Some  modifications  are  made  in  order  to  apply  these  methods  to  the 
same  problem.  In  section  7.1  and  7.2,  outlines  of  the  methods  by  Wiest  and  Tal¬ 
bot  are  given,  and  modifications  made  in  each  of  them  are  identified. 

The  two  modified  methods  and  the  intensity  assignment  algorithms  will  be 
used  in  scheduling  some  projects  and  results  will  be  compared.  The  most  satis¬ 
fying  results  would  be  obtained  if  the  project  data  were  real.  But  we  have  two 
limitations  in  using  real  project  data  in  comparison.  First,  it  is  very  hard  to  find 
a  reasonable  number  of  real  projects  which  will  give  a  fair  number  of  statistics 
for  comparison.  Second,  large,  realistic-size  projects  are  too  big  to  be  handled 
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by  Talbot's  method  efficiently.  Consequently,  arbitrarily  generated  small  pro¬ 
jects  will  be  used  in  comparison  of  methods.  In  section  7.3  ,  the  explanation  of 
the  arbitrary  network  generation  is  given. 

Four  methods,  Upgrading-Only,  Upgrading  and  Downgrading,  Wiest’s  heuris¬ 
tic,  and  Talbot's  Implicit  Enumeration,  are  used  to  schedule  120  arbitrarily  gen¬ 
erated  projects.  Results  are  presented  in  section  7.3.  Comparison  of  results  are 
made.  Note  that,  the  methods  by  Viest  and  Talbot  are  not  originally  designed  to 
solve  the  problem  with  our  assumptions,  so  .  the  comparison  of  results  is  only 
aimed  to  provide  evidence  of  satisfactory  performances  of  the  intensity  assign¬ 
ment  algorithms,  not  to  find  the  dominating  method. 

7.1.  Review  of  Methods  Used  in  Comparison 

7.1.1.  Wiest's  Heuristic  Model 

Viest's  heuristic  model  focuses  on  available  resources  like  intensity  algo¬ 
rithms  do.  It  serially  allocates,  period  by  period,  to  activities  available  to  start. 
Activities  are  scheduled,  starting  with  the  first  period,  by  selecting  from  the  list 
of  those  currently  available  and  ordered  according  to  their  total  slack.  The 
activity  with  the  least  slack  has  the  highest  priority  of  being  scheduled  as  avail¬ 
able  resources  permit. (Note  that,  SPAR1  uses  a  probabilistic  dispatching  rule, 
in  which  the  least  slack  activity  has  the  highest  probability  of  being  scheduled. 
A  simple  priority  dispatching  rule  is  used  here  since  a  one  pass  solution  will  be 
used.  A  multi  pass  solution  was  originally  used  in  SPARl  to  cope  with  minimizing 
the  total  cost  of  the  schedule.)  If  an  available  activity  fails  to  be  scheduled  in 


that  period,  an  attempt  is  made  to  schedule  it  in  the  next  period.  Eventually, 
activities  postponed  become  critical  and  move  to  the  top  of  the  priority  list  of 
available  activities. 

The  basic  model  described  above  is  modified  by  a  number  of  additional 
scheduling  heuristics  which  are  designed  to  increase  the  use  of  available 
resources  and/or  to  decrease  the  length  of  the  schedule. 

Crtw  Size:  With  each  activity  is  associated  a  normal  crew  size,  ie,  the 
number  of  men  normally  assigned  to  the  activity;  a  maximum  crew  size,  or  the 
maximum  number  of  men  which  can  be  assigned  to  the  activity;  and  a  minimum 
crew  size,  or  the  smallest  number  of  men  which  can  be  assigned  to  the  activity. 
Note  that,  these  crew  sizes  correspond  to  unit  intensity,  intensity  upper  T>ouml. 
and  intensity  lower  bound  ,  respectively,  in  the  variable  intensity  model,  and 
they  will  be  used  accordingly  in  solutions  of  arbitrary  problems.  The  rules  for 
crew  size  selections  are  as  follows:  If  an  activity  to  be  scheduled  is  critical,  it  is 
placed  on  a  priority  list  and  given  special  treatment.  If  sufficient  men  are  avail¬ 
able,  the  activity  is  scheduled  at  its  maximum  crew  size.  If  insufficient  men  are 
available  to  do  so,  or  even  to  schedule  the  activity  at  normal  crew  size,  then  an 
attempt  is  made  to  obtain  the  required  men  by  means  of  'Borrow  and 
Reschedule'  routines  described  below.  If  all  efforts  fail,  and  the  activity  cannot 
be  scheduled  at  even  minimum  crew  size  then  early  start  date  is  delayed  one 
period.  Non-critical  activities  are  scheduled  at  normal  crew  size  if  the  required 
number  of  men  are  available,  but  if  the  resources  available  are  insufficient  for 


scheduling  even  the  minimum  crew  size  then  the  activity  is  delayed  for  con- 


sideration  until  the  next  period.  Note  that,  an  activity  is  considered  to  be  criti¬ 
cal  if  its  slack  does  not  exceed  a  value  of  k,  where  k  is  an  input  parameter.  Zero 
is  assumed  for  the  value  of  k  in  the  experiments. 

Augment  Critical  Activities:  Repeated  attempts  are  made  to  speed  up  criti¬ 
cal  activities  which  have  crew  size  less  than  their  maximum  crew  size.  Before 
any  new  activities  are  scheduled  on  a  given  period,  activities  previously 
scheduled  and  still  active  are  examined.  If  any  of  these  is  critical  and  has  a 
crew  size  assigned  less  than  its  maximum,  and  if  resources  are  available,  its 
crew  size  is  increased  up  to  the  maximum. 

Ifulti- Resource  Activities ■'  Separate  activities  are  created  for  each  resource 
and  activities  are  considered  to  start  at  the  same  period  and  with  the  same 
level  of  resource  assignment.  Note  that  the  variable  intensity  model  handles 
multi-resource  activities  without  creating  any  extra  activities. 

Barrow  From  Active  Activities:  If  resources  available  are  insufficient  for 
scheduling  some  critical  activity  j,  then  the  model  enters  into  a  procedure  for 
searching  currently  active  activities  to  see  if  sufficient  men  might  be  borrowed 
from  them  to  schedule  the  activity  j.  Men  are  borrowed  from  activities  only 
when  the  resultant  stretching  of  the  activities  will  not  delay  the  entire  project. 

Reschedule  Active  Activities:  If  sufficient  men  could  not  be  obtained  from 
the  borrow  routine  described  earlier,  then  the  model  scans  the  list  of  currently 
active  activities  and  picks  out  those  which  could  be  postponed  without  delaying 
the  project  finish  time.  If  sufficient  men  is  are  found  in  this  way,  then  activity  j 
is  scheduled  and  necessary  adjustments  are  made.  Step  3  of  the  intensity 


Assignment  algorithms  performs  a  similar  procedure. 

Add  -on  Unused  Resources:  After  as  many  activities  as  possible  are 
scheduled  on  a  given  day,  there  still  may  be  unused  resources  in  some  shops 
for  that  day.  The  model  compiles  a  list  of  active  activities  to  -which  these 
resources  might  be  assigned  and  arranges  them  in  ascending  order  of  their 
total  slack.  Proceeding  down  the  list,  the  model  increases  the  crew  size  of  these 
activities  until  the  unused  resources  are  fully  allocated  or  else  the  list  is 
exhausted.  A  similar  procedure  is  performed  in  Step  4  of  intensity  assignment 
algorithms. 

Viest's  model  is  designed  to  find  the  schedule  with  minimum  cost  for  the 
overall  project.  The  model  tries  to  generate  many  different  schedules  and  each 
schedule  is  evaluated  by  a  cost  function.  After  a  number  of  iterations  the 
minimum  cost  schedule  is  chosen  by  a  search  routine.  Since  our  interest  is  in 
minimizing  the  project  duration,  search  routines  will  not  be  used  in  experi¬ 
ments.  A  one-pass  solution  will  be  used. 

Viest's  model  basically  assumes  three  alternative  levels  of  resource  appli¬ 
cations  for  each  activity  (Normal  crew  size,  maximum  crew  size,  minimum  crew 
size).  In  his  article,  Viest  mentions  that  some  intermediate  crew  sizes  may  be 
selected  although  it  is  not  shown  how  it  is  done  in  the  flow  chart  presented  in 
the  article  (Viest  (1967)).  He  proposes  the  intermediate  crew  sizes  to  still  have 
discrete  values  for  the  number  of  men  per  day.  Our  version  of  Viest’s  model 
uses  the  variable  intensity  model  which  allows  a  continuous  range  of  resource 
applications  within  upper  and  lower  bounds.  Our  version  of  Viest's  model 


assumes  that  intermediate  levels  of  resource  applications  (crew  size)  can  be 
chosen  only  when  augmenting  the  critical  activities. 


Wiest’s  model  and  the  intensity  assignment  algorithms  have  similar  heuris¬ 
tics  for  exploiting  the  use  of  available  resources  and  decreasing  the  project 
length.  The  orders  of  application  of  those  heuristics  are  different  in  each 
method.  For  example,  West's  model  first  augments  all  the  critical  activities  in 
the  active  list  and  then  starts  scheduling  the  ready  activities.  The  intensity 
assignment  algorithm  puts  both  ready  and  active  activities  into  a  combined  list 
which  is  sorted  by  their  slacks  and  augments  from  this  list.  Earlier  experiments 
with  the  intensity  assignment  algorithms  suggest  that  using  a  combined  list 
produces  better  schedules  then  augmenting  first  the  active  list  and  then  the 
ready  list. 

A  FORTRAN  code  was  created  for  Wiest's  model.  Originally,  Wiest  used 
periodic  review  for  resource  assignments,  but  event-based  review  is  used  in  this 
code  which  is  more  efficient.  Resource  assignments  are  made  only  when  there  is 
a  change  in  the  system.  Continuous  resource  application  is  assumed. 

7.12.  Talbot's  Implicit  Enumeration  Method 

Talbot  (19B2)  discusses  a  model  which  is  more  general  than  the  classical 
resource  constrained  project  scheduling  problem.  The  model  permits  each 
activity  to  be  accomplished  in  one  of  several  different  modes.  Each  mode  has 
predefined  resource  requirements  and  duration.  A  schedule  specifies  how  each 
activity  should  be  accomplished  (in  which  mode),  that  is,  which  resource- 
duration  option  should  be  selected,  and  when  each  activity  should  start.  A  two 


phase  implicit  enumeration  methodology  is  proposed  to  solve  the  problem.  The 
procedure  proposed  gives  optimal  solutions  for  small  size  problems  but  works 
as  a  good  heuristic  for  larger  problems. 

Talbot's  model  assumes  resources  consumed  and  durations  of  activities  are 
discrete.  Once  a  mode  is  selected  for  an  activity,  resources  used  by  the  activity 
are  constantly  applied  from  its  start  to  its  finish.  This  is  quite  different  than  the 
assumptions  of  the  variable  intensity  model,  namely,  that  intensity  of  an 
activity  may  vary  between  its  bounds  during  the  execution  of  the  activity. 

In  order  to  apply  Talbot's  model  to  the  same  problem  with  intensity  assign¬ 
ment  algorithms,  some  assumptions  are  made.  It  is  assumed  that  all  modes  of 
an  activity  have  the  same  resource  requirement  mix  but  application  rates  of 
each  mode  are  different.  Each  activity  has  a  number  of  modes  of  operation  (Aj. 
The  number  of  modes  of  activity  i  is  calculated  as  follows: 

*  "  0.1/d*  1 

Intensity  upper  bound  of  the  activity  corresponds  to  fastest  mode  (mode  1)  and 
intensity  lower  bound  of  the  activity  corresponds  slowest  mode  (mode  Kt)-  An 
intermediate  mode  k  represents  an  intensity  level  s^  in  the  variable  intensity 
model,  which  is: 

-(*-l)0.1/d*  (7.1) 

Note  that  this  modification  of  Talbot's  model  still  has  discrete  alternatives  for 

resource  applications.  So  the  optimal  solution  for  this  case  may  not  be  the 

optimal  solution  for  continuous  variable  intensity  case. 
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Solution  Methodology 

Talbot  developed  a  two-stage  methodology.  Stage  one  defines  the  problem 
as  a  compact  integer  programming  problem,  while  stage  two  searches  for  the 
optimal  solution  using  an  implicit  enumeration  scheme  that  systematically 
improves  generated  heuristic  solutions. 

Stage  One 

In  stage  one  the  network  is  first  relabeled  using  one  of  the  scheduling 
rules.  Talbot  presents  results  for  eight  different  scheduling  rules  in  his  paper, 
and  the  Minimum  Late  Finish  rule  time  gave  the  best  performance  for  the 
minimum  project  duration  problem.  The  fastest  mode  for  all  activities  is 
assumed  in  the  late  finish  time  calculations.  In  our  experiments,  the  Minimum 
Late  Finish  rule  is  used  to  relabel  the  network.  The  effect  of  this  relabeling  pro¬ 
cedure  is  to  specify  the  order  in  which  activities  will  be  considered  for  assign¬ 
ment  during  stage  two  of  the  algorithm.  Resources  and  modes  are  also  sorted  in 
stage  one.  Resources  are  sorted  such  that  the  resource  having  the  maximum 
frequency  of  highest  per  period  requirement  relative  to  average  availability  has 
the  smallest  numerical  label.  The  effect  of  this  sort  during  the  enumeration 
procedure  is  to  identify  resource  infeasibility  as  early  as  possible.  Modes  are 
sorted  by  increasing  duration.  In  the  experiments,  each  activity  A%  has  K% 
modes  in  which  mt  corresponds  to  the  fastest  mode  (*»,(*)  =  *i). 
corresponds  to  the  slowest  mode  s  J«)  writh  the  rest  of  the  modes 

arranged  according  to  equation  7.1.  Note  that  the  sorting  of  the  modes  affects 
the  initial  project  completion  time  found  by  the  enumeration  procedure. 
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Stage  Two 

In  this  stage  an  implicit  enumeration  algorithm  builds  feasible  partial  solu¬ 
tions  into  complete  schedules  by  considering  activities  for  assignment  in 
increasing  numeric  order.  For  each  activity  modes  are  also  considered  in 
increasing  numerical  order,  ie,  fastest  mode  to  slowest  mode.  The  algorithm 
begins  with  an  attempt  to  find  the  earliest  feasible  assignment  of  activities  with 
labels  1,  2,  3.  and  so  on,  in  order,  until  an  activity  is  identified  that  cannot  be 
assigned  even  with  slowest  mode  without  violating  precedence  or  resource  con¬ 
straints  or  without  exceeding  a  bound  such  as  its  late  finish  time  LFj.  (A  large 
value,  ie.  sum  of  the  durations  of  all  activities  with  their  slowest  mode,  is  used 
as  the  late  finish  time  of  the  last  activity  A#  for  the  first  trial.  The  late  finish 
bounds  of  other  activities  are  calculated  by  the  CPM  method  backward  pass 
procedure  whereby  the  duration  of  each  activity  corresponds  to  its  fastest 
mode.)  When  such  an  activity  is  identified  the  algorithm  backtracks  to  the  most 
recently  assigned  activity  and  attempts  to  reschedule  it  at  a  later  time  in  order 
to  free  resources  for  use  by  the  unassigned  activity.  This  process  is  systemati¬ 
cally  applied  to  all  activities  and  modes  until  an  improved  solution  is  found  or 
until  optimality  is  verified.  If  an  improved  solution  is  found,  bounds  (  LFj's)  are 
correspondingly  tightened  (ie,  the  amount  of  difference  between  LF#  of  current 
trial  and  the  value  of  improved  project  duration  is  subtracted  from  the  current 
late  finish  bounds  of  activities  in  order  to  find  the  late  finish  bounds  of  activi¬ 
ties  that  will  be  used  at  the  next  trial.)  and  the  assignment  procedure  begins 
again  with  activity  1.  Ultimately  optimality  (or  termination)  is  verified  in  one  of 
two  ways.  Either  all  possible  activities  assignments  have  been  explicitly  or 


implicitly  evaluated  or  a  solution  is  found  which  equals  a  proven  theoretical 
bound  such  as  the  project  duration  with  minimum  possible  durations  for  all 
activities. 

A  FORTRAN  code  of  Talbot's  method  was  created.  In  our  experiments,  the 
number  of  trials  to  find  the  optimal  solution  is  restricted  to  five  trials.  The  com¬ 
puter  time  for  each  problem  is  also  restricted  to  60  seconds  of  CPU  time.  If 
optimality  could  not  be  reached  by  then,  the  existing  solution  (ie,  best  solution 
so  far)  is  taken  as  the  schedule. 

7.2.  Generating  Arbitrary  Project  Networks 

Data  for  a  reasonable  number  of  project  networks  for  comparison  of  pro- 

* 

ject  scheduling  methods  is  hard  to  find.  Project  networks  appear  in  such  a  wide 
range  of  different  contexts  that  no  one  class  of  problems  can  claim  to  be 
representative.  Thus,  it  was  decided  to  generate  and  solve  test  problems.  Test 
problems  were  generated  by  a  computer  model  which  generates  arbitrary  pro¬ 
ject  networks.  Project  networks  were  randomly  drawn  with  the  following  attri¬ 
butes: 

1.  Number  Of  Activities(Length) 

Number  of  activities  was  randomly  chosen  between  30  and  50.  This  size  of 
project  network  is  not  representative  of  the  size  of  projects  that  intensity 
assignment  algorithms  are  designed  for.  Talbot’s  method  can  solve  up  to 
this  size  of  problem  within  a  reasonable  computing  time. 

2.  Width  Of  Network 

This  is  the  maximum  number  of  arcs  found  in  any  vertical  cut  in  the 


network.  Random  manipulation  of  the  incidence  matrix  representing  pre¬ 
cedence  relations  of  activities  is  used  to  randomize  the  width  of  the  net¬ 
work.  The  detailed  description  of  the  procedure  is  given  in  Appendix  3. 

3.  Complexity  Of  Network 

This  is  the  ratio  of  the  total  number  of  arcs  to  the  total  number  of  nodes  in 
the  network.  It  is  randomized  between  1.5  to  4.5.  The  detailed  description 
of  the  procedure  is  given  in  Appendix  3. 

4.  Activity  Durations 

Activity  durations  are  randomized  between  0  and  10  units  of  time.  The 
range  is  chosen  small  because  the  efficiency  of  Talbot’s  method  depends  on 
the  time  horizon  chosen. 

5.  Resource  Use  Of  Activities 

It  is  assumed  that  there  are  five  types  of  non-storable  resources.  Each 
activity  can  use  up  to  three  different  types  of  resources.  The  types, 
amounts  and  number  of  different  resources  used  are  randomized  for  each 
activity. 

6.  Tightness  Of  Resources 

Resource  capacities  were  assumed  to  be  constant  through  each  project 
execution  in  all  experiments  for  simplicity.  Note  that  the  intensity  assign¬ 
ment  algorithms  can  handle  varying  resource  capacities.  SPARl  and  TIE 
methods  can  also  be  modified  to  handle  varying  resource  capacities  but 
the  examples  in  the  literature  for  both  of  these  methods  use  only  constant 
resource  capacities.  A  percentage  of  the  peak  capacity  level  of  each 
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resource  in  the  Early  Start  schedule  of  a  project  is  used  as  the  capacity  of 
that  resource.  Peak  capacity  levels  of  resources  used  in  the  Early  Start 
schedule  of  the  project  are  calculated  under  the  assumptions  of  fixed  nor¬ 
mal  intensities  for  all  activities  and  the  resources  of  all  types  are  unlim¬ 
ited,  ie.  as  it  is  in  CPM.  Each  resource  capacity  is  found  by  multiplying  its 
peak  utilization  level  by  a  random  draw  from  the  [  BOX. 80%  ]  interval. 

7.  Resource  Application  Flexibility 

The  project  network  generated  up  to  this  step  is  used  to  generate  four 
parametric  problems.  Resource  application  flexibility  is  represented  by  a 
parameter.  Flexibility  ranges  [lower  intensity  bound,  upper  intensity 
bound]  are  assumed  to  be  equal  for  all  activities  in  each  experiment.  Each 
project  network  is  to  be  scheduled  under  the  following  flexibility  conditions 
for  activity  intensities. 


A. 

JU  =  1-0/  di 

,  *;  a  1.0/  dx 

for  all  i  (Inflexible) 

B. 

JU  =  0.8/ d* 

,  =  1.2/  d* 

for  all  i 

C. 

=  0.5/  dj 

.  s'  =  1.5/ d* 

for  all  i 

D. 

=  0.2/ dt 

.  *i  =  1.8/ dt 

for  all  i  (Most  flexible) 

7.3.  Experiments  with  Methods  and  Analysis  of  Results. 

Thirty  arbitrary  project  networks  were  generated  by  the  procedure 
explained  in  section  7.2.  Characteristics  of  the  arbitrary  projects  are  summar¬ 
ized  in  Table  7.1.  First,  projects  were  scheduled  by  four  methods  (  1. 
Upgrading-Only,  2.  Upgrading  and  Downgrading,  3.  SPARl,  4.  Talbot's  Implicit 


Problem 

Number 

:  Number 

! 

j  Activities 

Number 

of 

Resources 

Network  Resource 

Complexity  Tightness 

Critical 

Path 

Duration 

1 

36 

5 

2.54 

65 

B5 

2 

37 

5 

2.03 

64 

61 

3 

48 

5 

1.63 

77 

50 

4 

47 

5 

1.72 

64 

44 

5 

31 

5 

3.45 

64 

B1 

6 

44 

5 

1.91 

64 

119 

7 

41 

5 

1.57 

64 

71 

B 

37 

5 

2.11 

64 

B2 

9 

44 

5 

2.73 

60 

90 

10 

38 

5 

2.55 

62 

139 

11 

44 

5 

3.02 

69 

90 

12 

36 

5 

2.57 

69 

44 

13 

48 

5 

2.29 

61 

77 

14 

44 

5 

1.36 

69 

51 

15 

44 

5 

3.00 

69 

77 

16 

33 

5 

2.73 

63 

98 

17 

39 

5 

2.3B 

70 

101 

18 

45 

5 

3.16 

77 

84 

19 

46 

5 

3.72 

76 

100 

20 

42 

5 

4.10 

74 

107 

21 

45 

5 

3.5B 

77 

64  ! 

22 

40 

5 

3.01 

78 

77 

23 

38 

5 

3.11 

62 

100 

24 

37 

5 

2.19 

74 

70 

25 

41 

5 

1.49 

66 

54 

26 

47 

5 

2.2B 

73 

79 

27 

37 

5 

2.46 

74 

110 

28 

30 

5 

3.63 

79 

105 

29 

36 

5 

1.94 

66 

49 

30 

36 

5 

2.14 

66 

51 

Average 

40.3 

5 

2.55 

68.7 

81 

Range 

30  -  4B 

5 

1.36  -  4.10  60  -  79 

44  -  139  | 

TABLE  7.1.  Project  Characteristics  of  artificially  generated  projects  which  are 
used  in  the  experiments. 


Enumeration  (TIE))  for  fixed  intensity  case,  ie,  such  that  classical  resource  con¬ 


strained  project  scheduling  problems  are  faced.  The  first  three  methods  pro¬ 
duced  the  same  schedule  for  each  of  the  projects,  since  they  use  the  same 
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priority  dispatch  rule  for  activity  starts  (Least  Slack  First).  The  TIE  method  pro¬ 
duced  different  schedules  for  some  of  the  projects.  Table  7.2  summarizes  the 
comparative  results  of  the  first  three  methods  against  the  TIE  method. 


Cases 

Frequency  Percentage 

Number  of  Projects  in  which 

All  Methods 

Produced  Equal  Schedules 

17  57 

Number  of  Projects  that 

TIE  Method  Produced 
Shorter  Schedule  than  Others 

7  23 

Number  of  Projects  that 

TIE  Method  Produced 
Longer  Schedule  than  Others 

6 

Total 

30  100 

Table  7.2  Summary  of  Project  Completion  Time  Results  For  Fixed  Intensity  Case 

[**1.0/4.4=1.0/41 

For  the  fixed  intensity  case,  the  TIE  method  finds  the  optimal  schedule 
unless  it  is  interrupted.  Interruption  may  happen  for  one  of  two  reasons:  either 
computer  time  allowed  for  the  problem  is  exceeded  or  the  number  of  trials 
required  to  shorten  the  schedule  exceeded  five.  In  the  experiments,  computer 
time  was  the  cause  of  interruption  for  all  of  the  problems  interrupted  before 
finding  the  optimal  schedule  for  this  method.  Results  show  that  in  23%  of  the 
projects,  the  TIE  method  could  not  find  the  optimal  solution  before  being  inter¬ 
rupted.  This  reflects  the  inefficiency  of  the  TIE  method  even  in  smaller  size 
problems.  In  the  fixed  intensity  case  there  is  only  one  mode  for  each  activity  so. 
the  problem  size  is  smaller  than  the  problems  with  variable  intensity  case. 

The  first  three  methods  produced  equally  good  or  better  schedules  than 


the  TIE  method  in  77%  of  the  projects.  This  reflects  the  fact  that  the  Minimum 


Slack  First  priority  rule  gives  good  results  for  the  classical  resource  con¬ 
strained  project  scheduling  problem.  This  supports  the  findings  of  Devis  and 
Paterson  (1975)  and  many  others. 

Next,  the  projects  in  Table  1  were  scheduled  by  the  four  methods  for  three 
resource  application  flexibility  cases  in  which  each  \  has  intensity  bounds  of 
[fcsO.e/d*  .  [fcsO.S/d*  .£-*  =  1.5/^]  and  [*=0.2/^  .  £-*  =  1.6/**]. 

respectively.  Project  completion  time  results  for  these  cases  are  reported  in 
Table  3.  Table  4.  and  Table  5,  respectively.  The  tables  include  the  percentages 
above  the  best  project  finish  time  by  each  method  for  each  project.  The  best 
finish  time  of  a  project  is  the  smallest  of  the  project  finish  times  found  by  the 
four  methods.  The  TIE  method  found  the  optimal  solution  for  some  projects. 
However,  recall  that  the  optimal  solution  found  by  the  TIE  method  may  not  be 
optimal  for  the  continuously  variable  intensity  case.  The  TIE  method  only  con¬ 
siders  a  number  of  discrete  values  of  intensities  for  the  activities,  and  cannot 
assign  any  value  between  upper  and  lower  intensity  bounds.  But  for  those 
cases  that  the  TIE  method  found  an  "optimal"  solution,  we  expect  the  solution 
would  be  near-optimal  for  the  continuous  intensity  case,  since  the  number  of 
discrete  intensities  that  the  TIE  algorithm  may  pick  was  chosen  reasonably 
large  (ie,  JQ  =  5  for  the  mild  flexibility  case,  Ai  =  ll  for  the  more  flexible  case,  and 
=  16  for  the  most  flexible  case)-t 

Table  7.6  provides  a  summary  of  project  completion  time  results  for  the 

three  cases  of  resource  application  flexibility.  Overall  results  show  that  the 

t  Thet  is,  the  TIE  method  would  be  near-optimal  for  the  continuous  intens.ty  cese  when  in¬ 
tensity  must  be  held  constant  from  start  to  finish  of  each  activity. 
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Upgrading  tc  Downgrading  algorithm  performed  best  by  giving  the  best  project 
completion  time  result  51  times  out  of  90,  averaging  3%  above  the  best  project 
completion  time,  and  producing  a  maximum  of  16%  above  the  best  completion 
time.  The  TIE  method  performed  second  best,  by  giving  the  best  or  equally  best 
results  47  times,  averaging  4.7%  above  best,  and  producing  a  maximum  of  27% 
above  best.  Third  and  fourth  best  are  Upgrading-Only  and  SPARl,  respectively. 

For  the  less  flexible  case  (  £4  =  0.8/  d*  ,  Fi  =  1.2/di  )  the  TIE  method  per¬ 
formed  a  little  better  than  the  Upgrading  and  Downgrading  algorithm.  For  this 
case,  the  number  of  modes  in  TIE  is  only  five  and  the  method  could  find  an 
optimum  more  frequently  than  in  the  cases  with  more  flexibility,  since  less  com¬ 
puter  time  is  required. 

For  higher  flexibility  cases  [£4  =  0.2/ =  1.0/ d,]  and 

Ui=0  .5 /  d%  ,  fi  =  1.5/d,].  the  Upgrading  and  Downgrading  algorithm  performed 
better  than  the  TIE  method  in  all  measures.  Comparing  Tables  7.3,  7.4  and  7.5 
indicates  that  increases  in  resource  application  flexibility  improves  the  perfor¬ 
mance  of  the  Upgrading  and  Downgrading  algorithm. 

The  Upgrading-Only  Algorithm's  performance  was  fairly  good  considering 
that  there  is  a  restriction  of  management  such  that  once  resources  are 
assigned  to  an  activity  they  can  not  be  withdrawn.  In  general,  the  Upgrading- 
Only  Algorithm  performed  better  than  SPARl.  Both  Upgrading-Only  and  SPARl 
become  noticeably  inferior  to  the  other  two  methods  when  one  increases  the 

i 

flexibility  of  resource  application.  Both  of  these  methods  are  more  prone  to  fall 
into  the  case  that  resources  are  tied  up  unnecessarily  by  some  nonurgent 

t 


Average  Percentage  I  Number  of  Problems  I  Maximum 
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activities  than  is  the  Upgrading  and  Downgrading  algorithm. 

Table  7.7  summarizes  the  relative  performance  of  the  methods  with  respect 
to  project  finish  time  results  of  the  120  experiments  (including  the  inflexible 
resource  application  case).  Results  again  indicate  that  the  Upgrading  and 
Downgrading  Algorithm  performed  the  best  of  all,  with  74  shortest  durations,  40 
of  which  were  the  unique  shortest  duration. 


METHODS 

Cases 

Upgrading 

Only 

Upgrading  8c 
Downgrading 

BUS 

Number  of  Times  Shortest 
Schedules  Produced 
Including  Ties 

28 

74 

24  72 

Number  of  Times 
Unique  Shortest  Schedule 
is  Produced 

0 

40 

0  46 

Number  of  Times 
Longest  Schedule 
is  Produced 

3 

0 

62  15 

i 

Table  7.7  Relative  Performance  of  methods  to  Project  Completion  Times  results 
for  all  experiments  including  inflexible  resource  application  case. 

In  none  of  the  test  problems  did  the  algorithm  provide  the  unique  longest  dura¬ 
tion  for  the  project.  Talbot's  method  gave  the  shortest  duration  72  times,  46  of 
which  were  the  unique  shortest  duration,  and  the  method  produced  the  unique 
longest  duration  18  times.  The  Upgrading-Only  Algorithm  found  the  shortest 
durations  28  times  including  ties.  In  none  of  the  cases  was  the  shortest  dura¬ 


tion  unique,  and  in  only  three  cases  did  it  produce  the  unique  longest  project 
duration.  With  respect  to  the  number  of  longest  durations  found  by  the 
methods,  Upgrading-Only  performed  better  then  TIE  method  by  3  to  18. 


7.4.  Computational  Experience 


Projects  were  scheduled  by  each  of  the  four  methods  with  an  IBM  3091 
computer.  Each  project  network  in  Table  7.1  was  solved  for  the  four  different 
resource  application  flexibilities  at  once  in  order  to  reduce  the  cost  of  the 
whole  experiment.  The  recorded  computation  times  were  for  the  batch  of  four 
problems  by  each  method.  Table  7.6  summarizes  the  computing  time  results  for 
the  batches  of  four  problems.  Talbot's  method  finished  the  batch  of  four  prob¬ 
lems  within  the  allowed  computer  time  16  times  out  of  30.  Allowed  time  for  a 
batch  of  four  problems  was  chosen  to  be  200  seconds  of  CPU  time.  In  the  cases 
that  the  200  second  time  limit  was  exceeded,  individual  runs  were  made  with  a 
60  second  time  limit  of  those  problems  which  were  not  finished  in  the  batch.  In 
this  way.  a  time  limit  of  60  seconds  for  each  problem  was  allowed.  The  results 
reported  below  for  Talbot’s  method  exclude  the  cases  that  the  program 
exceeded  the  allowed  time  limits. 


METh 

I0DS 

Cases 

Upgrading 

Only 

Upgrading  & 
Downgrading 

SPAR1 

1 

TIE  • 

Range  of 
Computing  Times 
(in  seconds) 

1.31  -  2.57 

2.09  -  6.91 

1.20-2.11 

0.45-  134.B 

Average  of 
Computing  Time 
(in  seconds) 

2.29 

4.27 

1.55 

26.25 

B3 


The  first  three  methods  have  stable  ranges  of  computing  times  and  reason¬ 
able  average  computing  times.  SPARl  has  the  lowest  average  of  1.55  seconds.  Its 
highest  computing  time  was  only  2.11  seconds  for  a  batch  of  four  problems. 
Upgrading-Only  and  Upgrading  and  Downgrading  have  average  computing  times 
of  2.29  and  4.27,  and  maximum  computing  times  of  2.56  and  B.91  for  a  batch  of 
four  problems,  respectively.  These  results  indicate  that  none  of  the  first  three 
methods  would  be  seriously  limited  by  computing  time  when  scheduling  larger 
projects,  ie,  projects  with  several  thousand  activities  and  hundreds  of 
resources. 

On  the  contrary,  the  computing  time  for  Talbot’s  method  varies  drastically 
for  different  problems.  Even  though  the  project  networks  were  chosen  to  be  of 
small  size,  the  program  could  not  finish  twelve  out  of  thirty  cases.  Excluding 
these  cases,  the  average  computing  time  was  26.25  seconds  and  the  maximum 
computing  time  was  134.6  seconds  for  the  batch  of  four  problems.  Experience 
indicates  that  computing  time  in  this  method  is  affected  by  problem  size,  net¬ 
work  complexity  and  resource  application  flexibility  (number  of  modes  for  each 
activity). 

The  computing  time  results  suggests  that  Talbot's  method  is  too  inefficient 
to  solve  (in  a  practical  manner)  problems  much  larger  than  the  problem  size  we 
chose  in  experiments.  Yet,  the  size  of  many  real  projects  is  several  multiples  of 
the  size  of  the  problems  picked  in  the  experiments. 


fiyjyj 


7.5.  Conclusion 


The  analysis  of  results  has  shown  that  the  Upgrading  and  Downgrading 
method  produced  better  results  for  project  completion  times  compared  to  the 
others.  Results  for  the  Upgrading- Only  Algorithm  are  fairly  good  considering 
that  an  additional  restriction  is  applied  by  management.  Results  also  indicate 
that  both  algorithms  are  computationally  efficient  for  scheduling  large  projects. 
The  TIE  method  produces  near  optimal  schedules  unless  it  is  interrupted.  Com¬ 
putational  experience  suggests  that  Talbot's  method  is  too  inefficient  to  handle 
problems  of  larger  size  real  project  networks. 

Both  SPARl  and  TIE  methods  were  not  originally  designed  for  handle  vari¬ 
able  intensities  for  the  activities.  These  methods  were  modified  for  the  pur¬ 
poses  of  comparison  of  results  of  these  two  methods  to  the  results  of  intensity 
assignment  algorithms.  The  validity  of  the  results  that  can  be  obtained  from 


these  algorithms  is  evident. 


B.  HANDLING  STOCHASTIC  ELEMENTS  IN  PROJECT  NETWORKS 


6. 1 .  The  Approach  for  Handling  Stochastic  Elements 

There  are  two  basic  types  of  uncertainty  in  project-oriented  production 
networks.  The  first  type  is  the  uncertainty  as  to  the  actual  amount  of  man¬ 
hours  required  to  complete  a  production  activity.  The  second  type  is  the  uncer¬ 
tainty  as  to  the  nature  and  amount  of  repair  rework  that  may  be  needed  after  a 
test  activity,  ie,  one  is  faced  with  a  stochastic  project  network. 

We  shall  use  Monte  Carlo  simulation  to  treat  the  stochastic  elements.  The 
reason  for  simulating  the  project  schedules  is  because  stochastic  elements  in 
the  project  cause  some  of  the  important  parameters  of  the  project  schedules  to 
be  random  variables,  such  as,  important  milestone  dates,  project  completion 
date,  and  total  use  of  resources  in  certain  periods.  The  probability  distribu¬ 
tions  of  such  random  variables  cannot  be  analytically  calculated  from  parame¬ 
ters  for  the  activities  of  the  project,  because  of  the  complex  convolutions 
involved  in  any  large  network.  By  repeated  simulation  of  randomly  selected 
values  for  the  stochastic  elements  in  the  network,  approximate  probability  dis¬ 
tributions  for  important  schedule  parameters  can  be  developed. 

The  methods  developed  to  handle  stochastic  elements  for  project  schedul¬ 
ing  in  the  literature  do  not  provide  treatments  directly  applicable  to  the  sto¬ 
chastic  elements  discussed  in  chapter  4.  Consequently,  a  simulation  model  is 
developed  which  handles  stochastic  resource  requirements  for  activities,  han¬ 
dles  stochastic  network  structure  of  the  project  and  uses  the  intensity  assign¬ 
ment  algorithms  to  exploit  the  flexibility  of  resource  application  to  the 


activities  which  is  available. 

The  approach  requires  modeling  the  structure  of  each  stochastic  element. 
Models  for  resource  requirements  of  production  activities  and  test  activities  will 
be  presented.  In  the  approach,  a  number  of  different  schedules  are  generated 
by  using  random  values  for  the  parameters  of  stochastic  elements  in  each 
simulation.  One  of  the  intensity  assignment  algorithms  will  be  used  as  a  decision 
maker  for  resource  allocation  among  activities.  The  resulting  schedules  from 
simulations  are  used  to  developed  approximate,  unbiased  distributions  for 
important  milestones  of  the  project  and  for  total  use  of  different  types  of 
resources  during  the  project  execution. 

B.2.  Stochastic  Resource  Requirements 

We  suppose  there  are  two  different  kinds  of  activities  classified  according 
to  their  stochastic  nature.  Those  in  the  first  type  are  called  "production  activi¬ 
ties”  which  include  manufacturing,  repairing  and  rebuilding  types  of  jobs.  Those 
in  the  second  type  are  called  "test  activities"  which  involve  the  performance  of 
tests  on  various  systems.  For  production  activities,  total  resource  require¬ 
ments  are  stochastic.  PERT  methods  attempt  to  model  this  by  assuming  that 
activity  durations  are  stochastic.  In  the  variable  intensity  model,  the  activity 
duration  is  a  function  or  resource  applications.  Stochastic  resource  require¬ 
ments  provides  a  better  representation  of  work  uncertainty  in  shipyards  than 
stochastic  activity  durations. 

Total  resources  required  by  an  activity  are  assumed  to  have  distributions 
with  known  parameters.  This  is  the  case  for  the  projects  of  a  shipyard.  The 
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shipyard  has  past  history  of  total  man-hours  from  each  resource  required  by 
activities.  Distributions  can  be  driven  from  the  histories. 

8.3.  Stochastic  Project  Network 

Resource  requirements  of  test  activities  are  fairly  stable  for  performing 
the  activity  itself,  but  every  test  activity  may  result  in  a  success  or  a  failure.  If 
it  is  a  success  then  its  regular  followers  may  start,  but  if  it  is  a  failure  ,  then 
there  needs  to  be  some  kind  of  repair  rework  utilizing  different  resources  than 
those  used  by  the  test.  The  repair  needed  may  change  with  the  type  of  failure 
in  the  particular  activity.  There  are  a  number  of  possible  repair  procedures 
(activity  subnetworks)  which  may  follow  each  test  activity.  After  one  of  the 
repair  procedures  is  completed,  the  test  activity  is  executed  again  with 
discounted  parameters  for  the  chance  of  further  rework.  Before  project  execu¬ 
tion,  it  is  uncertain  whether  a  test  activity  will  result  in  a  success  or  a  failure 
(or  even  repated  failures).  If  it  is  failure,  the  repair  subnetwork  to  be  under¬ 
taken  is  uncertain,  so,  the  network  of  the  project  is  stochastic. 

Model  For  Test  Activities 

Assume  At  is  a  test  activity.  Once  its  predecessors  are  finished  Ax  may  start 
and  perform  the  test  using  resources  (a*,).  The  result  of  Ax  can  be  either  suc¬ 
cess  or  failure.  If  it  is  success  its  regular  (deterministic)  followers  may  startf, 
but  if  it  is  a  failure,  one  of  the  repair  subnetworks  (/^)  must  be  undertaken. 
After  a  repair  subnetwork  is  completed  (A,)  is  to  be  performed  again.  A  network 

t  Strict  precedence  is  assumed  for  the  followers  of  test  activities. 
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representation  of  a  test  Ax  and  follow-on  rework  is  given  in  Figure  B.l.  In  addi¬ 
tion  to  previously  defined  variables,  test  activities  have  the  following  parame¬ 
ters; 

px  probability  that  test  At  is  a  success  at  the  first  trial. 

at  increase  in  probability  of  success  after  each  failure  of  test  A*.  (0<o*<l) 

Ti  discounting  factor  for  resource  requirements  after  each  failure  of  test  A*. 
(0<7t<l  ) 

*  ={*■••’■* . r«.} 

set  of  possible  repair  routes  tor  test  Note  that,  a  route  is  a  subnetwork 
of  activities.  For  each  activity  in  any  route,  resource  requirements,  dura¬ 
tion,  and  precedence  relations  must  be  specified.  Generally,  activities  in 
repair  routes  require  different  mixes  of  resources  than  the  test  activity 
requires. 

dru  normal  duration  of  repair  subnetwork  r4  after  test  Ax. 

q«  probability  that  repair  subnetwork  r4  is  executed  given  a  failure  of  test  Ax- 

£  *«  =  1.0 

t 


The  probability  that  test  Ai  is  success  (p4)  is  increased  by  a  constant  fac¬ 
tor  (l  +  Oj)  after  each  failure.  The  resource  requirements  of  are  discounted 
by  a  factor  (yx)  after  each  failure.  Given  a  failure  of  test  Ax,  the  probability 
that  a  repair  subnetwork  is  to  be  executed  (gw).  is  assumed  to  be  trial- 


invariant. 


0.4.  Ihitawnti  for  Sfcochaatic  Bmmto  la  rhnailatloaa 

At  the  beginning  of  each  simulation  of  the  project,  resource  requirements 
of  all  production  activities  are  randomized.  That  is, 

■m  *-  Oh  random* 

where  random^  is  a  random  number  drawn  from  a  known  multiplicative  distri¬ 
bution  for  Ai  with  mean  1..  Then,  randomized  resource  requirements  are  used 
as  actualized  resource  requirements  for  the  simulation. 

Test  activities  are  treated  differently  at  their  finish  times.  After  each  test 
the  following  procedure  is  performed. 

1.  Determine  if  the  test  is  a  success  or  a  failure,  by  a  random  draw. 

2.  If  a  failure,  determine  the  repair  route  to  be  taken,  by  a  random  draw. 

3.  Do  updates  for  next  trial  of  test  A, . 

a.  Discount  resource  requirements: 

«**  *■  7i  for  all  k 

b.  Increase  the  probability  of  success: 

Pi  Pi  +  t*i 

c.  Update  following  parameters  for  the  test  Ax. 

Ft(t)  *-  T  where  T  is  a  big  number  distinguishing  that  test  Ai  must  be 
repeated. 

»<{f)  =  0.  and  pcvt{t)  =  0. 

Since  the  probability  of  success  improves  after  each  failure,  in  a  finite 
number  of  trials  test  i  will  be  a  success  and  its  followers  may  start. 
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Considering  the  durations  of  rework  activities  which  follow  the  test 
activity,  the  overall  duration  of  testing  and  rework  can  fluctuate  greatly.  This 
raises  the  question  of  what  values  for  test  durations  to  be  used  in  (target)  late 
finish  time  calculations  which  are  used  by  intensity  assignment  algorithms.  We 
have  two  choices  for  test  durations.  One  is  ,  assuming  that  tests  will  not  fail, 
take  the  normal  duration  of  a  test  as  its  duration.  Second  is,  find  expected 
values  using  its  parameters.  The  first  approach  can  be  criticized  to  be  myopic. 
Since  there  is  a  substantial  risk  that  the  overall  testing  duration  may  be  multi¬ 
plied,  we  might  choose  misleading  target  times  associated  with  the  wrong  criti¬ 
cal  path.  As  for  the  second  approach,  using  expected  durations  in  CPIf  calcula¬ 
tions  has  been  shown  to  produce  overly  optimistic  schedules  in  the  literature 
(for  fixed  intensity  case)  (Clark  (1962),  MacCrimmon  and  Ryavec  (1964)).  On  the 
other  hand,  allowing  variable  intensities  for  activities  gives  us  more  correction 
possibilities  than  fixed  intensity  models.  Expected  durations  for  test  activities 
(including  rework)  are  used  in  late  finish  time  calculations. 

0.4.1.  Expected  Durations  Tor  Test  Activities 

The  expected  duration  for  a  test  activity  i  is  calculated  as  follows: 

•*  =  «*»♦(  F[AU  -  1  )  (  d*  ♦  ard.) 

where. 

E[A^]  Expected  number  of  trials  until  success  for  test  A *.  After  every  failure 


the  value  of  this  term  is  updated,  since  probability  of  success 
increases.  Initially, 

^[^]sP*+2(l-p,)(pt4at)4S(l-f>t)(l-(p<+o*))(p*42ai)+....4 


where  I*  is  the  maximum  number  of  unsuccessful  trials  that  test  A* 
may  have.  ie.  (I»*[(l-p4)/aJ+l)t-  After  k  unsuccessful  tests, 
•£[AUs(P*+*«<)+2(l-(p<+k<w))(Pi+(* +  !)«*)+•••+ 

erd*  Expected  repair  duration  following  an  unsuccessful  test  A*.  (=  jJ  qu  dru) 

i«i 

Note  that  expected  durations  of  test  activities  are  recalculated  after  each 
failure  and  the  most  recent  values  are  used  as  activity  durations  in  step  2  and 
step  3  of  the  intensity  assignment  algorithms. 

8.42.  Expected  Late  Finish  Updates 

Recalculations  of  expected  late  finish  times  of  activities  are  made  before 
every  assignment  of  intensities  in  the  simulations.  The  procedure  described  for 
late  finish  updates  in  section  6.1.3  can  easily  be  modified  for  expected  late 
finish  updates.  We  need  to  change  the  step  1(b)  of  the  procedure  in  section 
6.1.3.  Step  1(b)  will  be  changed  to  following  expression. 

(b)  If  the  activity  is  in  progress,  then 
AF%(t-l)  -  t-l  +  edt  -pevi(f-l)  d* 

Note:  for  a  production  activity  i  •di=dt,  so  that  for  production  activities, 
AFi(t-l)  =  <-l  +  (  1  -pcvtO-l) )  d* 

We  also  need  to  substitute  edj  for  d*  is  step2  of  the  procedure  in  section 
8.1.3. 

Vee  ■umt  at  is  chosen  so  thst  L f  Is  integer. 
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9.  SCHEDULING  A  SHIP  OVERHAUL  PROJECT 

The  model  developed  for  project  scheduling  in  project-oriented  production 
systems  was  used  to  schedule  a  ship  overhaul  project.  Data  of  the  project  was 
provided  by  a  naval  shipyard.  The  description  of  the  project  is  given  in  section 

9.1.  Several  experimental  schedules  were  generated  by  using  intensity  assign¬ 
ment  algorithms  under  different  assumptions  discussed  in  section  9.2. 

First,  ‘deterministic*  schedules  were  generated  for  two  different  resource 
application  flexibility  conditions.  The  deterministic  schedules  assume  that  the 
resource  requirements  of  the  activities  are  deterministic,  and  also  that  the  pro¬ 
ject  network  is  deterministic  .  Then,  stochastic  simulations  were  made  nsing  the 
Upgrading-Only  algorithm  as  the  decision  maker  for  the  resource  applications 
among  activities  (Section  9.2).  Then,  the  results  of  experimental  schedules  are 
compared  and  their  implications  for  the  shipyard  management  are  discussed  in 
section  9.3. 

9.1.  Ship  Overhaul  Project 

9.1.1.  Project  Activities 

The  project  network  describing  the  ship  overhaul  consists  of  2931  activi¬ 
ties.  They  are  classified  as: 

(i)  Those  requiring  resources  with  deterministic  followers,  ie,  production 
activities  (  Type  1) 

(ii)  Those  requiring  resources  with  stochastic  followers,  ie,  test  activities  (Type 

2) 


(iii)  Those  requiring  no  resources  and  that  have  fixed  duration,  ie,  time  res¬ 
traints  (Type  3) 

(iv)  Dummy  activities  for  precedence  relations,  ie.  those  that  have  zero  dura¬ 
tions  (Type  4) 

The  number  of  activities  of  each  type  were: 

Type  1  202? 

Type  2  249 

Type  3  198 

Type  4  457 

Activities  of  Type  3  were  required  to  establish  the  consistency  between  the 
due  dates  specified  by  the  aggregate  planning  level  and  the  activities  of  the 
project  network.  The  number  of  activities  in  the  deterministic  network  were 
2465,  and  the  number  of  repair  rework  activities  were  466. 

The  average  number  of  successors  to  an  activity  was  2.36,  indicating  sub¬ 
stantia]  connectedness  of  the  activities  in  the  network.  The  number  of  overlap¬ 
ping  precedence  relations  (ie,  flow  transfers)  between  activities  was  536,  which 
was  eight  percent  of  the  total  precedence  relations. 

9.12.  Resources,  Activity  Durations.  Resource  Requirements 

The  service  of  38  different  shops  are  used  to  carry  out  the  ship  overhaul. 
Each  shop  has  a  number  of  work  centers  which  correspond  to  a  specific  type  of 
trade  labor,  ie,  a  nonstorable  resource  type.  Table  9.1  tabulates  a  list  of  impor¬ 
tant  shops.  The  total  number  of  shop-work  center  pairs,  ie,  the  number  of 
resource  types,  was  193.  Each  activity  may  require  services  from  a  number  of 
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SHOP  NO. 

SHOP  NAME 

13 

Quality  Assurance-  Engineering  and  Analysis 

33 

Quality  Assurance-  Non-nuclear  Inspection 

34 

Quality  Assurance-  Metallurgy  and  Material  Control 

33 

Quality  Assurance-  Non- destructive  Testing 

30 

Production  Department 

Structural  Group-  Shipfitters 

17 

Structural  Group-  Sheet  Metal 

26 

Structural  Group-  Welders 

41 

Structural  Group-  Boilermakers 

81 

Structural  Group-  Foundry 

06 

Mechanical  Group-  Tool  Shop 

23 

Mechanical  Group-  Forge 

31 

Mechanical  Group-  Inside  Machinist 

38 

Mechanical  Group-  Machinist  Marine 

56 

Mechanical  Group-  Pipefitter 

51 

Electrical  Group-  Electrical 

67 

Electrical  Group-  Electronics 

64 

Service  Group-  Woodmakers 

71 

Service  Group-  Painters 

72 

Service  Group-  Riggers 

99 

Temporary  Service 

32 

Nuclear  Engineering 

Table  9.1.  The  list  of  important  shops  that  the  overhaul  project  needs  service 
for  execution. 

resources.  The  maximum  number  of  resources  required  by  an  activity  was  70. 
and  the  average  number  of  different  resources  required  by  an  activity  was  9.7. 

An  activity  may  contain  a  number  of  elementary  work  packages.  For  each 
work  package  of  an  activity,  a  time  duration  is  estimated.  A  work  package  is 
called  a  Key-Operation.  The  estimated  duration  of  an  activity  is  the  maximum  of 
the  durations  of  its  Key-Operations.  In  the  case  of  Dummy  activities,  the 
estimated  duration  is  zero,  and  for  activities  of  Type  3.  activity  durations 
correspond  to  fixed  time  lags. 

Total  shop  service  hours  required  for  an  activity  are  estimated  by  summing 
over  the  shop  service  hours  required  by  the  Key-Operations  of  the  activity.  The 
shop  service  hours  for  activities  of  Type  3  and  Type  4  are  zero. 
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Tbe  size  of  the  ship  overhaul  project  may  be  appreciated  by  adding  the 
service  man-days  over  all  activities  of  each  shop  as  illustrated  in  Table  9.2. 

9.1.3.  Bounds  On  Activity  Intensities 

In  a  practical  case  when  data  is  available  concerning  work  space  and 
structure  being  worked  on,  minimum  and  maximum  shop  hours  per  unit  time 
can  be  estimated  directly  for  each  activity.  However,  for  this  study  the  data  on 
practical  limits  were  not  available,  and  bounds  on  the  activity  intensities  had  to 
be  determined  from  suppositions. 

A  uniform  upper  bound  of  1.2/  if*  was  taken  for  activity  intensities,  ie.  a 
daily  allocation  of  shop  hours  twenty  percent  greater  than  needed  to  complete 
the  work  in  the  estimated  "normal"  time  (<0  would  not  be  exceeded.  A  uniform 
lower  bound  of  0.5/  rf*  was  taken  for  the  activity  intensities,  ie,  an  activity  can 
not  be  extended  more  than  200  percent  of  its  original  duration.  In  this  way  a 
daily  allocation  of  shop  hours  can  not  be  less  than  fifty  percent  of  the  shop 
hours  needed  to  complete  the  work  in  tbe  estimated  "normal"  time.  These 
bounds  apply  only  to  Type  1  and  Type  2  activities,  as  Type  3  and  Type  4  activi¬ 
ties  have  fixed  intensities,  ie,  upper  and  lower  bounds  are  equal  to  1  /d*. 

9.1.4.  Limitations  on  Resources 

The  aggregate  planning  level  provided  man-hour  availabilities  of  each  type 
of  resource  (Shop-Work  center),  ie,  preallocated  shop  capacities  for  the  project. 
Han-hour  availabilities  are  given  as  tbe  daily  man-hour  availabilities  from  each 
shop-work  center  pair.  The  daily  availabilities  change  at  the  beginning  of  each 
month  (*'•20  working  days). 

9.1.5.  Milestones 

Tbe  events  that  are  important  in  project  executions  are  called  mile'stones. 


Ten  milestones  of  interest  in  the  ship  overhaul  are  as  follows: 

(1)  31B  Start  steaming  Plant  2 

(2)  321  Start  steaming  Plant  1 

(3)  825  Start  Cold  Operations  in  Plant  2 

(4)  82B  Complete  Cold  Operations  in  Plant  2 

(5)  833  Complete  Hot  Operations  in  Plant  2 

(8)  635  Start  Cold  Operations  in  Plant  1 

(7)  838  Complete  Cold  Operations  in  Plant  1 

(8)  843  Complete  Hot  Operations  in  Plant  1 

(9)  918  Production  Work  Complete 

(10)  999  All  Services  Complete 

9.2.  Deterministic  Experiments 

Four  deterministic  schedules  were  generated  for  the  conditions  given 
below.  Deterministic  schedules  were  generated  by  assuming  the  resource 
requirements  of  activities  are  certain  and  equal  to  the  expected  man  hour  ser- 
vice  requirement  from  each  work  center,  and  by  assuming  the  test  activities  will 
be  successes  at  their  first  trials,  ie,  no  rework  activities  are  actualized. 
Schedules  were  generated  under  the  following  conditions: 

(1)  Resources  are  unconstrained  ,  ie,  unlimited  resource  capacity.  Resource 
applications  are  inflexible,  ie,  activities  have  fixed  durations  ( g,  =  £,  =  1  /  et, ) . 
Note  that  CPU  assumes  these  conditions;  the  schedule  generated  here  is 
equivalent  to  the  Early  Start  Schedule  in  the  CPM  method. 

(2)  Resources  are  constrained  according  to  aggregate  plans,  and  resource 
applications  are  inflexible  for  activities.  Note  that  the  classical  resource 
constrained  problem  assumes  these  conditions  (although  most  solution 


99 


procedures  for  the  classical  problem  cannot  cope  with  time-varying  capa¬ 
cities).  Both  intensity  assignment  algorithms  yield  the  fame  schedule  in 
this  case  since  they  use  same  priority  rule. 

(3)  Resources  are  constrained  and  resource  applications  are  flexible,  ie, 
Ms  =  0.5/  di  and  ff<  =  1.2/ rf*.  For  these  conditions  the  project  is  scheduled 
differently  by  the  two  intensity  assignment  algorithms;  Upgrading-Only 
Algorithm  (3A)  and  Upgrading  and  Downgrading  Algorithm  (3B). 

Results  of  actualization  of  important  milestone  dates  are  summarized  in 
Table  9.3 


Milestones 

Unconstrained 

Fixed 

Intensity 

(1) 

Constrained 

Fixed 

Intensity 

(2) 

Constrained 

Variable 

Intensity 

Upgrading 

Only 

(3A) 

Constrained 
Variable 
Intensity 
Upgrading  ft 
Downgrading 
(3B) 

31B 

265 

455 

315 

293 

321 

265 

486 

322 

305 

825 

277 

450 

352 

334 

828 

319 

480 

377 

360 

833 

355 

355 

404 

398 

835 

291 

450 

326 

329 

838 

335 

475 

350 

352 

843 

371 

371 

388 

393 

918 

305 

530 

331 

332 

999 

525 

553 

531 

529  1 

Table  9.3.  The  summary  of  results  for  milestone  actualization  times  for  deter¬ 
ministic  experiments. 

Comparison  of  schedule  1  and  schedule  2  indicates  that  there  are  some  tight 
resources  in  schedule  2  which  cause  the  delay  in  actualizations  of  all  mile¬ 
stones  with  respect  to  actualization  dates  obtained  from  schedule  1  in  which 
resources  are  unconstrained. 


Comparison  of  schedule  2  with  schedule  3A  and  3B  indicates  that  flexibility 
in  resource  applications  improves  the  schedule  for  the  milestone  actualization 
dates  generated  by  both  intensity  assignment  algorithms.  When  resources  can 
be  applied  in  variable  rates,  scheduling  under  the  assumption  of  fixed  durations 
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and  fixed  resource  application  rales  inhibits  productivity.  Comparison  of 
schedule  3A  with  schedule  3B  indicates  that  having  a  managerial  restriction 
over  resource  applications,  ie.  resources  can  not  be  withdrawn  from  an  activity 
once  they  are  committed  until  the  activity  finishes,  also  causes  increases  in 
milestone  dates. 

0.3.  Stochastic  Simulations 

To  perform  stachastic  analysis,  total  man  hour  requirements  of  each 
activity  from  each  work  center  was  assumed  to  be  a  random  draw  from  a  uni* 
form  distribution  between  BOX  and  120X  of  its  estimated  value,  since  the  distri¬ 
butions  based  on  historical  man  hour  requirements  were  not  available  for  this 
study. 

Repair  rework  subnetworks  following  249  test  activities  were  developed  by 
the  shipyard.  Estimated  probabilities,  resource  requirements  ,  and  durations 
for  repair  activities  were  also  provided. 

The  Upgrading-Only  intensity  assignment  algorithm  was  used  as  a  decision 
maker  for  resource  allocation  between  overhaul  activities  in  all  simulations. 
The  Upgrading-Only  algorithm  was  chosen  because  current  decision  making  for 
allocation  of  resources  between  activities  in  the  shipyard  can  be  best 
represented  by  this  policy. 

40  different  schedules  were  generated  by  simulations.  At  each  simulation, 
man  hour  service  requirements  were  randomized,  and  a  different  subnetwork 
after  each  test  activity  was  actualized  due  to  randomizations.  The  statistics 
calculated  considering  all  simulations  give  distributions  for  some  important 
project  parameters.  Table  9.4  gives  the  cumulative  probability  distributions  of 
actualizations  of  milestone  dates.  Table  9.5  provides  distributions  of  man  hour 
usage  histories  through  project  life  for  an  example  shop  (Shop  3B).  Table  9.6 
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provides  the  distributions  for  total  man-hours  consumed  by  ell  project  activi¬ 
ties. 

0.4.  Conclusion 

The  deterministic  experiments  have  shown  that  using  flexibility  in  resource 
applications  improves  the  schedule  for  milestone  actualization  dates  generated 
by  both  algorithms.  The  shipyard  actually  applies  flexibility  in  resource  applica¬ 
tions  to  the  project  activities.  Methods  that  assume  fixed  resource  application 
rates  for  activities  may  give  overly  conservative  schedules  for  overhaul  projects 
of  the  shipyard.  Thus,  shipyard  management  should  apply  a  schedule  which  is 
generated  by  a  method  that  assumes  flexible  resource  applications  for  project 
activities  as  its  guidelines  for  scheduling  the  activities  and  projecting  the 
schedule  parameters. 

The  deterministic  experiments  have  also  shown  that  increase  in  flexibility 
in  resource  applications  to  the  project  activities  improves  the  schedule  for 
milestone  actualization  dates.  In  actual  practice,  the  shipyard  management 
generally  does  not  allow  withdrawal  of  resources  from  an  activity  before  it  is 
finished  although  it  allows  supplementing  the  resources  to  speed  up  the  activity 
execution.  This  is  the  policy  which  is  used  in  Upgrading-Only  Algorithm.  This 
policy  has  advantages  for  the  shop  managers,  ie,  it  avoids  complexity  concern¬ 
ing  labor  reassignments.  The  disadvantage  of  this  policy  is  that  it  may  result  in 
longer  durations  for  projects.  The  shipyard  management  should  consider  the 
resource  assignment  policy  in  order  to  decide  if  the  policy  is  beneficial  to  the 
shipyard.  The  benefit  of  finishing  the  projects  earlier  may  be  worth  tolerating 
more  managerial  difficulties  at  the  shop  level.  The  schedule  generated  by  the 
Upgrading  and  Downgrading  algorithm  has  shown  that  the  policy  change  wil) 
reduce  the  milestone  actualization  dates. 
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Table  9.4.  Cumulative  Probability  distributions  of  milestone  actualization  times. 
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Stochastic  simulations  have  been  used  to  obtain  unbiased  distributions  for 
important  project  statistics.  The  simulation  model,  by  using  flexible  resource 
application  rates  in  deciding  the  resource  allocations  among  activities,  by 
characterizing  the  stochastic  man-hour  requirements  for  production  activities, 
and  by  characterizing  the  stochastic  network  structure  of  test  activities,  gave 
'good*  estimates  for  distributions  for  project  schedule  parameters.  Having 
'good'  distributions  for  project  parameters  will  improve  the  project  planning 
and  scheduling  process  in  all  levels  (inter-project  planning,  project  planning, 
project  control)  in  the  shipyard. 
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APPENDIX  A1  Restriction  on  Intensity  by  e  Progress  L*g 

According  In  the  definition  of  transfer  coefficient  / jt  ,  et  every  intensity 
assignment,  activity  intensities  must  be  assigned  so  that  the  progress  lags  from 
predecessor  activities  are  not  violated. 
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Figure  A.l  Intermediate  Product  Transfer 


Figure  A.l  presents  two  overlapping  activities,  A*  and  Aj  ,  Aj  is  a  predeces¬ 
sor  of  Ax  .  For  the  variable  intensity  case  equation  (5.5)  must  be  satisfied  for 
the  flow  transfers.  Rewriting  the  equation  (5.5)  after  substituting  cumulative 
intensity  terms  (pcvt)  gives 


P»,(0< 


Wax 

I 

1.0 


[o'pci/jCO  -/vj  if  OepnijU )  <1.0 
if  pcvj(t)=1.0 


For  the  case  pcv;(f)=  1.0,  choosing  intensity  for  Ax  is  independent  from  its 
predecessor  Ay  For  the  case  of  0*pcvj(t  )<1.0  the  following  inequality  must 
hold. 

pcv4 (t )  *  Max  |o  .  pci,(t  )~/jx 


The  algorithms  do  not  consider  assigning  intensity  to  Ax  before  time  ft, 
where  f,  is  such  that  pcvt(t  We  need  to  consider  only  the  following  case. 


pcVi(t)  *  pcvj(t)  -ffi  t  or 

Suppose  that  the  intensity  *x(t )  assigned  to  A*  will  be  maintained  until  t  +1. 
Then  the  following  inequality  must  hold: 

pcvt(t)  +  *t(f)  *pcvt(t)  -Ji i  +  *j(0 
Rearranging  this  gives 

**(0  ifc pcvj(t) -pcVi(t)  -ffi  +*j(t)  (Al.l) 

This  is  the  actual  restriction  on  intensity  of  A*  at  time  t.  This  assumes  that 

the  intensity  of  A*  will  be  changed  at  the  beginning  of  each  time  period.  We  do 
not  want  to  change  the  intensity  assignments  every  period  for  only  flow 
transfer  reasons.  Instead,  we  will  assign  a  conservative  intensity  level  to  Ax  such 
that  it  will  satisfy  the  progress  lag  requirements  for  all  periods  assuming  that 
the  assigned  intensity  of  A*  and  current  intensity  of  A;  will  be  kept  constant 
until  Aj  finishes.  Derivation  of  the  conservative  (more  restricted)  limit  on 
intensities  by  flow  transfers  is  as  follows. 

For  the  unit  intensity  case  ,  when  the  following  inequality  holds,  then  ine¬ 
qualities  (5.4)  and  (5.5)  both  hold. 

/;(t) * 

Generalizing  this  inequality  for  the  varying  intensity  case,  assume  that  activity  j 
will  keep  its  intensity  level  at  t  until  it  finishes  (  Fj(t )  ),  and  also  assume  the 
intensity  level  decided  for  activity  i  will  be  kept  constant  until  the  activity 
finishes  (  F«(f )  ),  Then  we  require 

r,(t) -MO*/*  1^77 


Since, 


10B 


1. 


by  substituting, 


Fi(t)  =  t  +  (1  -pcvi(t)) 


,4(1  -WO) 


Rearranging, 


(A1.2) 


.(,)**•=  P0-WO-/*) 

- FjiD-t 

This  must  hold  for  all  predecessors  Aj  of  A». 

The  equation  (Al.2)  imposes  a  tighter  restriction  on  intensity  of  A%  than 
the  equation  (Al.l)  does,  as  can  be  seen  from  the  following  argument. 

Since 


by  substituting  this  expression  into  the  equation  (Al.2)  and  rearranging  gives 


. -*«,(<)  -  /„)  *,(<> 

*,(,)* - r.o-Pwfi) 

From  (Al.l)  and  (Al.2),  we  have  to  show  the  following  inequality  holds 


pcvj  (f  )-~pcVi  (t )  +*>(f  )fc 


(1.0  -pCVj{t)  -fti)  Kj(t) 
1.0  -pcvj(tj 


Rearranging, 


(A1.3) 


(pcvj(t  )-pcVi(t)-/it)(l.O  -pcvjU))*  (pc*j(l)-pcvt(t)  -Sj i)  «j(0 
Rearranging, 


(1.0  -  pcvj(t  ))  2  Zj(t) 

This  inequality  holds,  because  the  intensity  of  Aj  .  ie,  fraction  of  Aj  which  will  be 
covered  in  period  t.  cannot  be  greater  than  the  uncovered  part  of  the  activity 
Aj-  This  is  true  by  definition.  Then  the  equation  (A1.3)  holds. 


r 


9 


Figure  A.2  demonstrates  the  difference  between  the  restrictions  imposed 
by  the  equation  (A1.2)  and  (A1.2). 


AC  represents  the  progress  curve  if  a*  is  chosen  as  intensity  level  (  Eq 


(A1 .2)).A§  ,  PC  represents  the  actual  limitation  imposed  by  the  progress  lag 
definition  (  Eq  (Al.l)). 
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APPENDIX  A Z  Handling  Storable  Resource* 

Suppose  that  some  of  the  resources  are  limited  in  their  total  (cumulative) 
consumption  as  well  as  their  usage  per  period.  Let 

Hi  ( t )  Consumption  limit  of  resource  type  1  up  to  time  t.  ie.  cumulative  con¬ 
sumption  of  resource  type  1  at  time  t  cannot  exceed  this  limit.  It  is  an 
increasing  function. 

[4  f  a*  *t(r)<  Aj(t)] 

T«1  t*l 

Then,  we  need  to  modify  the  intensity  assignment  algorithms  to  handle  this 
case.  Define: 

ARt(t)  Available  resource  amount  from  resource  type  1  at  time  t,  ie.  unused 
cumulative  resource  limit  at  time  t. 

*  ruo-  £  £ «« *»(■»■) 

<*1  r*l 

Basically,  the  feasibility  of  cumulative  resource  consumption  must  be 
satisfied.  This  can  be  done  by  checking  if  there  are  available  resources  to  con¬ 
sume  for  the  activities  in  the  next  period. 


The  following  modifications  in  intensity  assignment  algorithms  are  needed 
to  handle  this  type  of  resource  limitation: 

(1)  Initialize  the  cumulative  available  resources. 

AJ?i  =  J?,(l)  for  all  l 


(2).  In  step  2(c),  in  deciding  new  *t(f),  we  need  to  consider  the  availability  of 
cumulative  resource  consumptions.  The  equation  for  yt(f )  will  be  changed 


*  *•  y\  y*«  t*«  ,*•  *%  t*»  y'  kN  y**  «,**>*»  .*•  »  *  ,  •  •  '  «  <  ^  *  w  ^  >4  'j*  '  *  *  a  •  »  *  »  *  *  *  i'*  »  ’  *  *  ,  *  »  •  \  '%**'**» 


TOW 


to  the  following  form. 


Vi(< )  =  Uin 


min 

* 


f(  M tensity  Upper  Bound, 

ark(t) 


min 


M*(0 


+  *i(t  -1)  Rssaurct  Availiblity 


i  | — a - Rssourcs  Consumption  Availablity 


(l.O-j»ei/t(f-l)) 


1.0 


Tbrgst  Finish  Tims 


mto  |(1.0-pcvi(t— 1)— /ji)  flow  Trims fwr 

then  update; 

ARl(t)  +  ARt(t)  -«,*(**(<  )-s*(<-l))  tor  all  1  if**(f)>0. 

(3)  In  step  3(c-i).  The  resources  for  the  rescheduled  activities  must  be  with* 
drawn. 


AR{  (t )  =  ARi  (t )  +  £  aj j  s*(t)  for  all  l  for  all  jtlRESCHED 


T*1 


(4)  In  step  4  ,  the  same  change  shown  tor  step  2(c). 

(5)  In  step  5d.  following  update  will  be  added 

A^(t)-A^(t-1)  -  a*  *«(<).  r=t . tn 

and  if  ARt( r)  becomes  negative  for  any  t,  then  apply  capacity  change  rou¬ 
tine  in  order  to  recover  the  infeasibility. 

(6)  Change  in  limit  of  consumption  will  be  treated  as  an  event  which  is  called  a 
*  consumption  limit  change  The  event  occurs  when  /fc(f)  >  R*(t- 1)  for 
any  k.  When  there  is  such  an  event,  the  following  updates  are  done.  (Let 
A/4(0  =  R*(t)  -/fc(f-D) 
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APPENDIX  A3  Generating  Arbitrary  Project  Networks  Using 
Randomized  Incidence  Matrices 

The  first  step  of  artificial  network  generation  was  the  generation  of  an 
arbitrary  incidence  matrix  that  represents  the  precedence  relations  of  project 
activities.  The  incidence  matrix  has  N  rows  and  N  columns,  where  N  is  the 
number  of  activities  in  the  project  which  is  randomized  between  30  and  50. 

Let  Xy  be  the  element  of  incidence  matrix  at  ith  row  and  jth  column. 

fl  if  At  is  followed  by  Aj 
“  jo  Otherwise 

First,  we  set  all  elements  of  the  incidence  matrix  to  1.  (Later  we  will  put 
some  elements  to  0  to  eliminate  the  precedence  relations  that  are  not  possible.) 
Then,  the  possible  relations  will  be  randomized  in  order  to  get  the  final  form  of 
the  incidence  matrix. 

The  project  networks  are  acyclic;  therefore,  the  part  of  the  matrix  that  is 
under  the  diagonal  must  be  set  to  zero. 

Xy  » 0  if 

In  project  networks,  it  is  unlikely  to  have  precedence  relations  between  the 
activities  which  are  at  the  beginning  of  the  network  and  the  activities  which  are 
at  the  very  end  of  the  network.  In  order  to  avoid  such  a  possibility,  one  can  put 
a  cut  which  is  parallel  to  and  above  the  diagonal  of  the  matrix.  The  column 
where  the  cut  starts  determines  the  width  of  the  network  (M).  The  column  M 
where  the  cut  starts  is  randomized  within  [(N/2)-5,  (N/2)+5].  Once  M  is  decided, 


the  cut  will  be  defined  as  follows, 
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Xij  s  0  if  U*j-i 

The  resulting  incidence  matrix  is  shown  in  Figure  A3.1.  The  elements 
within  the  region  ABCD  represent  the  possible  precedence  relations.  The  next 
step  will  be  to  randomly  define  zero  elements  within  this  region. 


D 


Figure  A3.1.  Project  network  incidence  matrix 

The  number  of  followers  (F)  for  each  activity  is  randomized  between  1.5  to 
4.5.  Note  that  most  of  the  project  networks  that  are  used  in  the  literature  fall 
within  this  range.  This  is  done  as  follows.  The  maximum  possible  number  of  fol¬ 
lowers  of  Ai  is  Ki  ,  where 

Ki*ZX* 

i 

In  order  for  A*  to  have  F  followers  on  average,  we  choose  the  following  probabil¬ 
ity  level  for  each  possible  follower  of  A*: 


Then,  a  random  draw  is  made  for  each  possible  follower  of  A^ie,  for  each 


entry  in  row  i  which  lies  in  region  ABCD).  Let  m  be  a  random  number  from  a 
uniform  distribution  within  [0,1].  Then. 


5 


Xy  *  0  if  rm>Pri 

This  procedure  is  performed  for  all  activities  of  the  project  in  order  to  find 
the  final  form  of  the  incidence  matrix  which  will  be  taken  as  precedence  rela- 
tions  of  the  arbitrary  project  networks.  Note  that,  the  expected  number  of  fol¬ 
lowers  F  is  not  the  same  as  the  network  complexity.  The  network  complexities 
that  are  reported  in  Table  7.1  are  calculated  after  incidence  matrixes  are  final- 
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